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SUMMARY 


Defects  In  semiconductors  lead  to  poor  reliability  and  loss  of  yield  of 
electronic  devices  fabricated  from  these  materials.  Currently,  x-ray 
topography  Is  used  to  characterize  semiconductor  wafers.  However,  the  long 
time  needed  to  characterize  each  sample  by  this  technique  precludes  its  use 
for  each  Individual  sample.  A  technique  thai  is  fast,  accurate,  efficient 
and  able  to  characterize  each  sample  in  a  reasonable  time  Is,  hence, 
desirable.  One  such  technique  utilizing  infrared  piezo-birefringence  is 
Investigated  in  this  work.  The  technique  is  based  on  changes  In  the 
refractive  Index  of  the  material  under  stress.  The  stress  could  be  from  an 
applied  external  load  or  due  to  defects  in  the  material  generated  during  its 
growth  and/or  processing.  The  purpose  of  this  research  is  to  Investigate  the 
use  of  infrared  birefringence  In  defect  characterization  of  semiconductor 
materials. 

As  a  first  step  towards  defect  image  simulation,  the  case  of 
diametrically  loaded  discs  of  semiconductor  materials  is  examined.  This  is 
done  to  obtain  a  better  understanding  of  the  simulation  algorithm  prior  to 
its  subsequent  use  In  defect  characterization.  A  dark-field  plane 
polarlscope  was  constructed  with  a  He-Ne  laser  tuned  to  1.15  pm  wavelength 
used  as  the  light  source.  This  particular  wavelength  was  selected  as  the 
semiconductor  materials  of  primary  Interest,  namely  SI  and  GaAs,  are 
transparent  at  this  wavelength  in  the  near  infrared.  The  computer  simulated 
Images  matched  remarkably  well  with  the  experimentally  observed  ones. 

The  behavior  of  the  stress-optic  coefficient  C  was  also  investigated, 
which  has  been  treated  as  a  constant  by  other  investigators  in  earlier 

works.  In  this  work.  It  was  found  that  for  (100)  oriented  Si  and  GaAs  discs 

under  diametrical  compression,  the  stress-optic  coefficient  C  is  a  strong 
function  of  position  for  a  given  load,  and  also  changes  with  the  orientation 
of  the  applied  load  with  respect  to  the  principal  crystal  axes.  However,  no 
such  dependence  was  found  for  (111)  oriented  SI  and  GaAs  discs  under 
diametrical  compression,  as  expected  from  the  crystal  symmetry.  The  values 

of  C  for  (100)  oriented  SI  disc  under  diametrical  compression  ranged  from 

2.0  x  10-12  cm2/dyne  to  3.0  x  lO-*2  cm2/dyne  and  for  (1001  oriented 
GaAs  disc  under  diametrical  compression  ranged  from  0.8  x  10~12  cm2/dyne 
to  2.6  x  10-12  cm2/dyne.  The  corresponding  figures  for  (111)  oriented  Si 
and  GaAs  discs  under  diametrical  compression  are  2.33  x  10~12  cm2/dyne 
and  1.94  x  10"12  cm2/dyne,  respectively. 

Next,  an  algorithm  was  developed  for  simulation  of  images  of  edge  and 
screw  dislocations  in  semiconductors.  For  simplification  purposes,  image 
simulations  were  carried  out  for  edge  and  screw  dislocations  when  the 
dislocation  line,  the  Burgers  vector  and  the  direction  of  the  incident  light 
were  along  the  crystal  principal  axes.  For  this  special  case  when  the 
incident  light  was  perpendicular  to  the  Burgers  vector,  the  computer 
simulation  found  distinguishing  features  between  an  edge  and  a  screw 
dislocation.  When  the  polarization  angle  of  the  incident  light  with  respect 
to  the  dislocation  slip  plane  was  45",  the  image  of  a  screw  dislocation  was 
found  to  disappear.  Whereas,  for  the  same  polarization  angle,  the  image  of 
an  edge  dislocation  was  brightest  with  a  four-petal  rosette  pattern  with 
equal  lobes. 


An  extension  of  this  work  Is  being  carried  out  In  our  laboratory  to 
generalize  the  algorithm  for  screw  and  edge  dislocations  to  take  Into  account 
any  arbitrary  orientation  of  the  dislocation  line  and  the  Burgers  vector 
alongwlth  any  arbitrary  direction  of  the  Incident  light.  In  principle,  a 
variety  of  these  Images  for  varying  polarization  angles  of  the  Incident  light 
can  be  generated  and  stored  on  a  computer.  A  sample  with  an  unknown  defect 
can  be  placed  under  the  polarl scope  and  the  Images  monitored  for  varying 
Incident  light  polarization  angles  and  Its  Incident  directions.  Then  an 
Image  matching  procedure  can  be  applied  to  match  the  experimentally  generated 
Images  with  the  simulated  ones  stored  In  the  computer.  A  direct  comparison 
will  then  give  the  type  of  defect.  Its  location,  and  Its  orientation  with 
respect  to  the  principal  crystal  planes. 


FOREWORD 


This  research  task  was  started  with  Dr.  Craig  S.  Hartley  as  the  project 
director  In  August,  1984.  However,  due  to  a  change  In  plans.  Dr.  Hartley 
left  LSU  In  August,  1986.  Since  then,  the  project  was  carried  out  under  the 
direction  of  Dr.  Pratul  K.  Ajmera. 
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Figures 


2.1.1  The  sample  geometry.  Axis  X3  Is  the  outward  normal  from  the 

plane  of  the  paper .  4 


2.2.1  Sample  with  x-z  plane  as  the  dislocation  slip  plane. 


A  right-handed  screw  dislocation  along  the  axis  of  a  cylinder 
of  radius  R  and  length  L . 


(a)  The  general  view  of  the  experimental  set-up.  (b)  The 
schematic  arrangement  of  the  experimental  set-up . 


(a)  Experimentally  observed  image  on  a  (111)  silicon  disc  of 
radius  0.873  cm  and  thickness  0.16  cm.  Diametrical  compres¬ 
sion  load  P  =  19  kg,  angle  of  polarization  a  =  90°.  (b)  Simu¬ 

lated  Image  for  the  same  conditions  as  (a).  Normalized  Inten¬ 
sity  Is  plotted  on  a  logarithmic  scale  for  the  computed  image.. 

(a)  Experimentally  observed  and  (b)  computed  Images  on  a  (111) 
silicon  disc  of  radius  0.873  cm,  thickness  0.16  cm  for  a  load 
of  38  kg  and  the  angle  of  polarization  a  of  90*.  Normalized 
Intensity  Is  plotted  on  a  logarithmic  scale  for  the  computed 
image . 

(a)  Experimentally  observed  Image  reported  by  Appel 
et.  al.  [11]  for  a  silicon  sample  of  unknown  orientation 
with  a  radius  of  0.437  cm,  thickness  of  0.188  cm,  load  of 
38  kg  and  the  angle  of  polarization  a  of  the  Incident  light 
of  60°,  (b)  simulated  Image  for  the  same  conditions  as  (a) 
for  (100)  orientation,  and  (c)  same  as  (b)  but  for  (111) 
orientation.  Normalized  Intensity  is  plotted  on  a  logarithmic 
scale  for  the  computed  Images . . 


(a)  Simulated  image  on  a  (100)  gallium  arsenide  disc  of 
radius  0.375  inch  and  thickness  65  mils.  Diametrical  com¬ 
pression  load  P  *  20  kg,  angle  of  polarization  a  *  0°. 

(b)  Experimentally  obtained  Image  for  the  same  conditions 

as  (a).  Normalized  intensity  Is  plotted  on  a  logarithmic 
scale  for  the  simulated  image . 


(a)  Simulated  image  on  a  (100)  gallium  arsenide  disc  of 
radius  0.375  inch  and  thickness  65  mils.  Diametrical  com¬ 
pression  load  P  =  40  kg,  angle  of  polarization  a  =  45°. 

(b)  Experimentally  obtained  Image  for  the  same  conditions 

as  (a).  Normalized  Intensity  Is  plotted  on  a  logarithmic 
scale  for  the  simulated  image . 
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3.4.1 


3.4.2 


6.1.1 


6.2.2 


6.2.3 


Wafer  map  of  the  computed  values  of  the  stress-optic  coeffi¬ 
cient  C  for  a)  S  -  45®  and  b)  S  =  75°  for  a  load  of  19  kg, 
(100)  silicon  wafer  with  radius  of  0.873  cm  and  thickness 
of  0.16  cm.  Plotted  in  11  levels  of  gray  scale  with  white 
for  C  >  2.95  x  10"12  cm2/dyne  and  black  for 
C  <  2.05  x  10"12  cm2/dyne . 


Wafer  map  of  the  computed  values  of  the  stress-optic  coeffi¬ 
cient  C  for  a)  5  =  45®  and  b)  6  =  75®  for  a  load  of  20  kg, 
(100)  gallium  arsenide  wafer  with  radius  of  0.375  Inch  and 
thickness  of  65  mils.  Plotted  In  10  levels  of  gray  scale 
with  white  for  C  >  2.5  x  10-12  cm2/dyne  and  black  for 
C  <  0.9  x  10-12  cm2/dyne . , 


Simulated  fringe  patterns  for  an  edge  dislocation  on  a  (100) 
Si  surface  with  the  dislocation  slip  plane  along  one  of  the 
principal  crystal  plane  and  the  Burger's  vector  along  one  of 
the  principal  crystal  direction,  a  Is  the  angle  between  the 
dislocation  slip  plane  and  the  polarization  axis  of  the  in¬ 
cident  light,  a)  a  =  0®,  b)  a  =  30®,  c)  a  =  45®,  d)  a  =  60®. 

Simulated  fringe  patterns  for  a  screw  dislocation  on  a  (100) 
SI  surface  with  the  dislocation  slip  plane  along  one  of  the 
principal  crystal  plane  and  the  Burger's  vector  along  one  of 
the  principal  crystal  direction,  a  Is  the  angle  between  the 
dislocation  slip  plane  and  the  polarization  axis  of  the  in¬ 
cident  light,  a)  a  =  0®,  b)  a  =  25®,  c)  a  =  50®,  d)  a  =  75®.. 

The  loading  and  the  crystal  principal  axes  systems  for  a 
(100)  plane . 


6.2.1  The  loading  axes  coordinate  system. 


The  crystal  axes  system  and  its  orientation  with  respect  to 
the  plane  of  the  wafer . 


The  (x^,  x£,  xp  coordinate  system  and  its  relationship  to  the 
loading  axes  system  and  the  crystal  principal  axes  system . 


88 


A 


m 


5*5' 


viii 


Tables 


LIST  OF  TABLES 


Page 


2.1.1  Conversion  technique  from  tensor  to  reduced  notation .  9 

2.1.2  The  piezo-optic  coefficients  for  Si  used  in  computation 

(after  [2]) .  9 

2.1.3  The  piezo-optic  coefficients  of  GaAs  as  derived  from  data  of 

Booyens  and  Basson  [9] .  9 

6.3.1  The  values  of  the  elasto-optic  coefficients  (Pij)  and  the 
stiffness  (c-j j )  for  GaAs  after  Booyens  and  Basson  [9]. 

(Pij  are  dimensionless,  c^ j  are  in  dyne/cm2.) .  51 


1.  INTRODUCTION 


Defects  In  semiconductor  crystals  play  a  major  role  In  the  final 
yield  of  a  manufacturing  process,  where  the  yield  Is  defined  as  the  number  of 
good  chips  on  a  wafer  divided  by  the  total  number  of  chips  on  the  wafer.  The 
defects  also  have  a  significant  impact  on  device  properties  and  parameters 
such  as  minority  carrier  lifetime,  carrier  mobility,  carrier  generation  rate. 
Impurity  diffusion,  optical  properties  -  such  as  the  Index  of  refraction  and 
the  optical  absorption  spectrum.  These  factors  can  lead  to  poor  reliability 
of  the  device  and  low  yield.  An  efficient  way  of  characterizing  these 
defects  will  not  only  Improve  the  yield  of  the  manufacturing  process,  but 
also  Improve  the  device  performance. 

The  most  commonly  used  technique  today  for  defect  characterization  Is 
X-Ray  topography  [1].  The  theory  behind  the  technique  Is  Bragg's  law  of 
diffraction,  generalized  by  Laue.  However,  the  relatively  long  exposure  time 
needed  for  characterizing  one  wafer  makes  the  turn-around  time  quite  long  (of 
the  order  of  1-2  days).  This  fact  precludes  testing  each  individual  sample, 
and  the  usability  of  a  lot  is  determined  by  statistical  estimates.  The  fact 
that  all  the  wafers  are  not  tested  allows  some  defective  wafers  to  pass 
through.  Therefore,  a  technique  which  is  fast,  accurate,  and  efficient,  and 
which  can  characterize  each  sample  in  a  reasonable  time  is  likely  to  improve 
the  yield  and  device  performance.  The  technique  being  considered  here  is 
Infrared  piezo-birefringence,  which  may  prove  useful  in  fast  and  accurate 
characterization  of  defects  in  semiconductor  materials. 

Most  semiconductor  crystals  such  as  silicon  and  gallium  arsenide 
employed  for  fabricating  electronic  devices  are  optically  Isotropic.  Both 
these  materials  are  opaque  to  visible  light  but  are  transparent  in  the  near 
infrared.  Silicon  having  a  band  gap  of  1.12  eV  has  an  absorption  cutoff  at 
1.11  pm.  The  corresponding  values  for  gallium  arsenide  are  1.42  eV  and 
0.89  pm,  respectively.  Light  with  wavelength  longer  than  1.11  pm  will 
normally  pass  through  these  two  materials.  The  light  does  not  experience  any 
relative  phase  change  while  traversing  through  these  optically  Isotropic 
samples.  Therefore,  if  a  sample  Is  placed  between  the  polarizer  and  the 
analyzer  with  their  axes  crossed,  a  dark  image  will  be  seen.  However,  under 
stress,  the  sample  becomes  blrefrlngent.  The  stress  can  be  applied 
externally  or  generated  Internally  due  to  defects,  processing,  or  material 
preparation.  The  light  while  passing  through  the  stressed  regions  of  the 
sample  splits  into  two  components  In  mutually  orthogonal  directions  having 
different  refractive  Indices  and  thus  velocities  along  the  two  directions.  A 
phase  difference  builds  up  between  the  two  components  as  they  traverse  the 
sample.  Now,  If  the  sample  Is  viewed  through  a  crossed  polarizer  and 
analyzer  set-up,  an  image  will  be  seen.  The  Image  will  be  characteristic  of 
the  phase  shift  between  the  two  components.  For  each  kind  of  defect,  the 
stress  fields  will  be  different  and  thus  different  fringe  patterns  or  Images 
will  be  seen  for  different  kinds  of  defects  such  as  an  edge  dislocation  or  a 
screw  dislocation.  This  phenomenon  of  stress-induced  splitting  of  the 
refractive  Index  Into  two  components  Is  known  as  piezo-birefringence. 

Glardlnl  [2]  studied  the  phenomenon  of  birefringence  In  silicon  and 
obtained  the  values  of  Its  piezo-optic  coefficients  relating  the  stresses  to 
the  changes  In  the  refractive  Indices.  Lederhandler  [3]  and  Hornstra  and 


Penning  [4]  observed  residual  stresses  In  silicon  by  the  Infrared 
birefringence  method.  These  stresses  were  attributed  to  severe  thermal 
gradients  present  during  crystal  growth  which  resulted  In  frozen-in  strain  In 
the  material.  Jenkins  et.al.  [5-6]  have  studied  Images  of  dislocations  In 
silicon  samples  by  the  piezo-birefringence  method.  The  stresses  produced  by 
the  dislocations  cause  the  normally  optically  Isotropic  silicon  to  become 
blrefrlngent. 

In  the  work  carried  out  here.  Infrared  piezo-birefringence  Is  utilized 
to  examine  Image  patterns  obtained  by  a  dark  field  plane  polarlscope 
from  Infrared  transparent  semiconductor  samples  under  stress.  The  normally 
optically  Isotropic  semiconductor  samples  become  blrefrlngent  due  to  the 
stress.  In  this  work.  Image  patterns  for  defects  have  been  simulated  on  a 
computer.  In  some  cases,  the  simulated  Images  are  compared  to  the 
experimentally  observed  ones.  As  a  first  step  to  obtain  defect  Images  and  to 
obtain  a  better  understanding  of  the  computer  Image  simulation  algorithm; 
Images  were  simulated  for  externally  applied  loads  on  SI  and  GaAs  discs.  The 
experimentally  observed  and  simulated  Images  were  compared.  Section  2 
describes  the  Image  simulation  algorithms.  Subsections  2. 1-2.3  discuss  the 
algorithms  for  semiconductor  discs  under  diametrical  compression,  edge 
dislocations,  and  screw  dislocations,  respectively. 

Section  3  outlines  the  experimental  work  done  and  discussions  about  the 
results.  Subsection  3.1  gives  the  experimental  set-up  utilized  In  this 
work.  Subsection  3.2  Illustrates  the  results  for  dlametrlcall.  compressed 
semiconductor  discs.  The  experimentally  obtained  Images  are  compared  with 
the  simulated  ones.  In  our  work,  we  Investigated  the  behavior  of  the 
stress-optic  coefficient  C,  which  was  taken  to  be  a  constant  In  earlier 
works  by  other  Investigators.  The  results  of  our  findings  are  given  In 
Subsection  3.3.  Subsection  3.4  discusses  the  simulation  results  for  edge  and 
screw  dislocations.  Section  4  gives  the  summary  of  our  work  along  with  the 
conclusions. 


2.  IMAGE  SIMULATIONS 


2.1  Semiconductor  Discs  under  Diametrical  Compression 


In  this  section.  Image  simulation  is  carried  out  for  Infrared 
transparent  semiconductor  discs  under  externally  applied  load  In  the  form  of 
diametrical  compression.  This  Is  done  to  obtain  a  better  understanding  of 
the  simulation  algorithm  and  as  a  first  step  towards  obtaining  simulation 
Images  for  defects  In  semiconductors.  ! 

i 

The  disc  under  consideration  is  of  radius  R,  thickness  d  and  under  a  I 

diametrical  compression  load  P.  The  loading  axes  system  Is  defined  In  j 

Fig. 2. 1.1.  The  loading  axis  X2  Is  chosen  arbitrarily  and  d  '.s  not  i 

necessarily  coincide  with  any  of  the  crystal  principal  axes.  Axis  xi  Is  ' 

perpendicular  to  axis  X2  and  lies  on  the  sample  surface  and  axis  X3  Is  I 

coming  out  from  the  plane  of  the  paper.  The  stress  at  any  point  (xi,X2)  ! 

for  a  disc  under  diametrical  compression  Is  given  by  [7] 
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2  2  2  2  2  2 
where  r^  =  x^  +  (R-x2)  and  r2  =  +  (R+x2)  .  Here  0^  is  the  stress 

produced  In  direction  x^  transmitted  across  a  surface  perpendicular  to 
direction  Xj.  The  on  and  022  components  are  the  normal  components  of 
stress  while  012  is  the  shear  component  with  012  =  021  because  there  is 
no  resultant  torque  In  the  system.  In  order  to  compute  the  stress  components 
Ojj,  the  sample  Is  divided  into  small  squares  and  the  stresses  are 
calculated  at  the  center  of  each  of  these  squares.  The  semiconductor 
materials  of  primary  Interest  In  this  work  are  Si  and  GaAs  both  of  which 
belong  to  the  cubic  class  of  crystal  symmetry.  The  Indlcatrlx  of  a  cubic 
crystal  under  unstressed  condition  Is  spherical  with  no  preferred  optic 
axis.  However,  under  stress  the  Indlcatrlx  changes  to  an  ellipsoid  of 
revolution  and  can  be  given  by  a  quadric  surface  [8] 
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where  k  and  1  are  summation  Indices  and  are  the  elements  of  the 

dielectric  Impermeability  matrix  defined  by  +  AB^,  where  B^  are 

the  elements  of  the  dielectric  Impermeability  tensor.  For  an  Isotropic 
material,  the  only  non-zero  values  of  are  given  by  B^  =  B22  =  B23  =  BQ, 


where  BQ  -  and  where  nQ  Is  the  isotropic  value  of  the  refractive  index. 
no 

Here,  ABy  are  the  changes  In  By  due  to  the  stress  components.  When  the 

light  is  incident  along  -X3  direction,  the  intersection  ellipse  as  obtained 
from  Eqn. (2.1.4)  Is  given  by: 

(Bo+ABll)xl2  +  (BQ+AB22)x2^  +  ^B12x1x2  +  &B21X2X1  “  1*  (2.1.5) 

The  diagonal Izatlon  of  Eqn. (2.1.5)  would  give  the  major  and  the  minor 
axes  of  the  ellipse  formed  by  the  intersection  of  the  Incident  light 
direction  with  the  Indlcatrlx.  The  result  of  the  diagonal Izatlon  yields  Bi 
and  B2,  the  dielectric  Impermeabilities  along  the  minor  and  the  major  axes 
of  the  Index  ellipse  respectively.  They  are  given  as: 


^Bii  +  1  ?  \l? 

B1  =  B0  +  — - -  +  5  C(AB1X-AB22)  +  4AB12AB21]  . 


(2.1.6) 


AB^  +  AB22 


\  C(AB11-AB22)2  +  4AB12AB21]] 


(2.1.7) 


The  corresponding  refractive  indices  are  given  by  nj  *  iNB\  and  ng  = 
I/VB2.  Due  to  these  local  differences  In  the  refractive  Indices,  the 
components  of  a  plane  polarized  light  beam  traversing  the  sample  will  show  a 
phase  shift.  When  the  light  emerging  from  the  sample  Is  viewed  through  an 
analyzer,  the  resultant  Image  will  show  a  pattern  characteristic  of  this 
phase  shift. 

The  piezo-optic  coefficients  are  used  to  compute  the  changes  In 

the  dielectric  impermeability  elements  AB^j,.  The  piezo-optic  coefficients 
are  conventionally  defined  along  the  crystal  principal  axes.  Therefore,  the 
computation  Is  carried  out  in  two  steps.  First,  the  stress  components  0^ 
obtained  from  Eqns. (2.1.1-2.1.3)  are  transformed  to  the  crystal  principal 
axes.  Then,  the  changes  In  the  dielectric  impermeability  elements  are 
obtained  In  the  crystal  principal  axes  system.  These  are  then  converted  back 
to  the  loading  axes  system. 

The  equations  for  transforming  the  stress  components  from  the  loading 
axes  system  (xltx2,x3)  to  the  crystal  principal  axes  system  (xj,x£,xp  are 

given  by  C8]  as 


ak1afcj°1 j 


(1-1, 2, 3;  j-1,2,3) 


(2.1.8) 


where  1  and  j  are  summation  Indices  and  a^  Is  the  direction  cosine  between 
the  crystal  principal  axis  x'|<  and  the  loading  axis  x^.  The  calculation 
of  the  direction  cosines  for  (100)  and  (111)  oriented  samples  under  arbitrary 
loading  direction  are  discussed  In  Appendices  6.1  and  6.2,  respectively. 
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In  terms  of  the  piezo-optic  coefficients  tiijkj,,  the  changes  In  the 
dielectric  Impermeability  AB'y  In  the  crystal  principal  axes  system  due  to 
the  loading  are  given  by  [8]  as 
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(2.1.9) 


where  k  and  1.  are  summation  Indices  and  o'^j,  are  given  by  Eqn.(2.1.8).  The 
piezo-optic  coefficients  uijkj,  form  a  fourth-rank  tensor  having  In  general 
81  terms.  A  reduced  notation  suggested  by  Nye  [8]  Is  often  used  to  express 
the  fourth-rank  tensor  coefficients  as  given  In  Table  2.1.1.  The  dielectric 
Impermeability  elements  AB'y  evaluated  by  Eqn.(2.1.9)  are  defined  along 
the  crystal  principal  axes.  The  last  step  of  the  procedure  Is  to  transform 
these  elements  back  to  the  loading  axes  system,  which  Is  given  by  the  trans¬ 
formation  [8] 
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where  1  and  j  are  summation  Indices  and  AB'^j  are  defined  by  Eqn.(2.1.9). 
The  Impermeability  elements  which  are  found  from  Eqn. (2.1.10)  are  used 
In  Eqns.(2.1.6)  and  (2.1.7)  to  obtain  the  major  and  the  minor  axes  of  the 
intersection  ellipse  formed  by  the  Incident  light  with  the  Indlcatrlx  sur¬ 
face.  Once  Bi  and  B2  are  known,  the  corresponding  refractive  Indices  can 
be  found  easily  from  nj  =  l/i/Bj  and  n2  -  I//B2,  where  nj  and  n2 
are  the  refractive  Indices  along  the  minor  and  the  major  axes  of  the  Index 
ellipse,  respectively.  After  the  refractive  Indices  are  known,  the  Intensity 
calculation  can  be  carried  out  as  given  below. 

A  light  wave  after  passing  through  the  polarizer  and  Incident  upon  the 
sample  has  the  form  A  cos(wt-$+2irx3/X)  where  X  is  the  wavelength,  w  Is  the 

angular  frequency  and  4  Is  an  arbitrary  phase  angle.  The  component  amplitude 
along  the  direction  of  the  minor  axis  of  the  Index  ellipse  Is  given  by 
A  cos(a-P).  As  shown  In  Fig. 2. 1.1,  a  and  p  are  respectively  the  angles  sub¬ 
tended  by  the  Incident  light  polarization  direction  and  the  minor  axis  of  the 
Index  ellipse  with  respect  to  the  axis  xj.  The  angle  p  Is  given  by  [7] 


P  =  |  tan”1 


AB^  +  AB22 
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Similarly,  the  component  amplitude  along  the  major  axis  direction  of  the  in¬ 
dex  ellipse  Is  given  by  A  sln(a-p).  The  light  Intensity  transmitted  through 
the  analyzer  will  be  polarized  along  the  analyzer  direction.  The  components 
along  the  major  and  the  minor  axis  directions  are  further  resolved  along  the 
polarization  and  the  analyzer  directions.  Let  the  component  amplitudes 
A  cos(a-P)  and  A  sln(a-P)  be  denoted  by  Aj  and  A2,  respectively.  Thus, 
for  a  ray  travelling  along  the  minor  axis,  the  component  amplitudes  along  the 
polarization  and  the  analyzer  directions  are  given  by  Aicos(a-P)  and 
Aisln(a-P),  respectively.  Similarly,  for  the  ray  travelling  along  the 
major  axis,  the  component  amplitudes  along  the  polarization  and  the  analyzer 
directions  are  given  by  A2Sln(a-P)  and  -  A2Cos(a-p),  respectively. 
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The  net  light  amplitude  L  coming  out  of  the  analyzer  Is,  hence,  written  as 
L  =  A  slnfa-PJcosfa-PJcosfwt-^ird/AjJ-A  s1n(a-P)cos(a-P)cos(wt-$+2*d/A2) 


where  A^  *  A/nj  and  A2  *  A/n2  are  the  wavelengths  along  the  minor  and  the 
major  axes  respectively.  The  above  expression  can  be  written  as 


L  =  (A/2)s1n2(a-0)[cos(wt-0+2ud/A1)  -  cos(wt-$+2*d/A2)]. 


Squaring  and  taking  the  time  average  to  get  the  net  resultant  Intensity  I 
coming  out  of  the  analyzer,  one  finally  obtains 


I  *  L2  =  (A2/2)  s1n22(a-p)  sln2{ird(l/A1-l/A2)}. 


(2.1.12) 


In  general,  for  each  point  on  the  wafer,  the  stresses  and  hence,  the 
angle  p  and  the  values  for  A*  and  A2  will  be  different.  Thus,  the 
resultant  intensity  image  will  show  dark  and  bright  regions  resulting  In  a 
fringe  pattern.  The  computer  program  developed  In  this  work  uses  the  above 
developed  analytical  approach  to  compute  the  Intensity  of  the  light 
transmitted  through  a  stressed  sample.  The  size  of  the  squares  for  the 
stress  calculation  was  selected  carefully  so  that  the  Intensity  did  not  show 
any  abrupt  discontinuity.  For  an  Idea  about  the  sizes  of  the  squares.  It 
will  suffice  to  say  that  to  simulate  the  Images  for  a  diametrically 
compressed  SI  disc  of  1.746  cm  diameter,  the  total  area  was  divided  Into 
approximately  60,000  squares.  Silicon  and  gallium  arsenide,  the  two  main 
semiconductors  of  Interest  In  this  work  belong  to  the  cubic  class  of  crystal 
symmetry;  silicon  being  a  member  of  m3n  group,  and  gallium  arsenide  being  a 


member  of  43m  group.  All  crystals  having  cubic  symmetry  have  only  three 
Independent,  non-equal,  and  non-zero  piezo-optic  coefficients  out  of  a 
possible  total  of  81.  In  reduced  notation,  these  three  non-zero  coefficients 
are  denoted  by  *n,  *12,  and  *44  [8].  Explicitly  they  can  be  expressed 
as  follows: 
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For  silicon,  the  values  used  In  Eqn. (2.1.9)  In  this  program  were  taken 
from  Glardlnl  [2]  and  are  shown  In  Table  2.1.2.  These  values  were  measured 
at  a  wavelength  of  1.11  pm  by  Glardlnl  [2]  which  Is  close  to  the  wavelength 
of  1.15  pm  used  In  this  work.  For  gallium  arsenide,  to  the  best  of  our 
knowledge,  the  values  for  the  piezo-optic  coefficients  are  not  available  in 
the  open  literature.  However,  the  paper  by  Booyens  and  Basson  [9]  has 
reported  the  elasto-optlc  coefficients  and  the  stiffness  values  of  gallium 
arsenide  at  a  wavelength  of  1.15  pm.  The  procedure  of  obtaining  the 
piezo-optic  coefficients  from  the  elasto-optlc  coefficients  and  the  stiffness 
values  Is  given  In  Appendix  6.3  and  the  calculated  values  are  shown  In 
Table  2.1.3.  In  the  computer  program,  quadruple  precision  arithmetic  Is  used 
throughout  to  maintain  numerical  accuracy. 

2.2  Edge  Dislocation 

This  section  deals  with  the  Intensity  distribution  around  a  straight 
edge  dislocation  viewed  by  Infrared  light.  The  sample  Is  taken  to  be  a 
square  with  the  slip  plane  perpendicular  to  the  sample  surface.  It  Is 
assumed  that  the  slip  plane  makes  an  angle  of  \|i  with  one  of  the  crystal 
principal  axes.  In  this  section,  as  before,  all  primed  notations  refer  to 
the  crystal  axes  system  and  all  unprimed  notations  refer  to  the  sample  axes 
system.  The  approach  given  here  Is  quite  general  and  is  demonstrated  here 
for  a  specific  situation.  As  a  specific  example,  we  take  a  (100)  surface  of 
silicon  crystal  with  the  dislocation  slip  plane  lying  perpendicular  to  the 
surface.  In  this  example,  we  assume  that  the  crystal  x'-axls  comes  out  of 
the  plane  of  the  paper  and  the  crystal  y'-axls  makes  an  angle  of  »1>  with  the 
dislocation  line  which  Is  also  taken  to  be  the  sample  x-axis  as  shown  in 
Fig. 2. 2.1.  This  assumption  Is  made  to  simplify  the  computations  for  this 
example. 

The  stresses  associated  with  the  dislocation  are  given  by  the  following 
formulae  [10]: 
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Table  2.1.2 

The  piezo-optic  coefficients  for  Si  used 
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Here,  p  Is  the  shear  modulus,  b  Is  the  magnitude  of  the  Burger's  vector  and  v 
Is  the  Poisson's  ratio.  Also,  021  =  012,  and  031  =  032  =  0,  as  there 
cannot  be  any  resultant  torque  In  the  system. 

Next,  the  direction  cosines  a^j  between  the  crystal  and  the  sample 
axes  systems  are  evaluated.  These  are  needed  to  transform  the  stress 
components  to  the  crystal  axes  system.  For  the  specific  example  chosen  here, 
they  are  given  by 

a^  =  cos  x^xj  *  cos  x'V  *  cos  90  =  0, 
aj2  =  cos  X2*x^  =  cos  yAx'  =  cos  90  =  0, 
a13  =  cos  x^xj  =  cos  z~x'  =  cos  0=1, 
a21  =  cos  x^x^  =  cos  x~y'  =  cos  \|i, 

822  =  cos  x2aX2  =  cos  y^y'  =  cos  (90-Y), 

823  *  cos  x3~x£  »  cos  z'V  =  cos  90  =  0, 
a31  =  cos  x^x^  *  cos  x^z'  =  cos  (90+x|») , 
a32  =  cos  x2ax^  =  cos  y"z'  «  cos  \|», 


a33  »  cos  x3~x3  x  cos  z~z*  =  cos  90  =  0. 


(2.2.6) 


Once  the  direction  cosines  are  evaluated,  the  stress  components  can  be 
transformed  to  the  crystal  axes  system  through  the  following  transformation: 


°kH  ak1a&j°1j  (i=l, 2, 3;  j=l,2,3), 


(2.2.7) 


Each  of  the  0^  term  Is  summed  over  nine  terms.  However,  023  ■  o32  »  o13  x 
o31  *  0,  as  well  as,  o12  *  o2r  Taking  these  Into  account,  Eqn. (2.2.7) 
simplifies  to 

°kl  '  aklaHl°ll  +  (aklaH2  +  ak2a!ll)012  +  ak2aH2°22  +  ak3afc3°33*  (2,2,8) 


k 

8 

9 

« 

& 

a 

a 


I 

s 

9 

v 

V 

V 


For  the  specific  case  under  consideration  here  for  which  aii=ai2=a23=a33=0, 
and  a.-  «  1,  Eqn.  (2.2.8)  can  be  given  by 


°il  =  °33* 

°12  *  °13  =  °21  =  °31  =  °* 

°22  *  a212°ll  +  2a21a22°12  +  *22°22' 

°23  "  °32  *  a21a31°ll  +  (a21a32  +  a22a31*°12  +  a32a22°22’ 
and 

°33  *  a312°ll  +  2a32a3I°12  +  a322°22* 


(2.2.9) 


The  changes  In  the  dielectric  Impermeability  tensor  elements  AB I ^  can  be  cal¬ 
culated  from  the  formula 


ABij  =  E  ^ijk^ka  (k=1‘2»3;  4=1,2, 3). 

k ,  Sr 


(2.2.10) 


The  cubic  crystal  symmetry  results  In  only  three  Independent  values  of  the 
piezo-optic  coefficients  vuw  as  given  by  Eqns. (2.1.13-2.1.15).  Taking 
this  Into  account,  the  coefficients  In  Eqn. (2.2. 10)  can  be  written  as 


ABil 

*  "llll0!!  +  vU22a22 

AB12 

"  B13  *  AB21  *  AB31  * 

AB22 

=  Tf2222°22  +  1r2211°il 

AB23 

=  AB'2  =  1T2323°23  +  *; 

and 

AB33 

"  w3333°33  +  *3311°il 

(2.2.11) 


These  dielectric  Impermeability  tensor  elements  are  defined  along  the  crystal 
axes.  These  are  transformed  back  to  the  sample  axes  through  the  following 
transformation 


(2.2.12) 


A8kjl  =  l  a1kajlAB'j  (1  =  1, 2, 3;  j=l,2,3).  (2.2.12) 

1 ,  j 

Each  of  the  AB  kj,  term  Is  summed  over  nine  terms.  Substituting 
Eqns. (2.2.11)  and  (2.2.6)  in  the  above  equation  yields 


ABn 

“  >212‘B22  + 

2a21a31AB23  + 

a31^B33’ 

AB^ 

=  AB2i  =  a2ia22AB22  +  a21a; 

32AB23  +  a, 

AB13 

“  AB31  °  AB23  =  AB32  =  °* 

ab22 

=  a222AB^2  + 

2a22a32AB23  + 

a32^AB33  * 

and 

ab33 

=  a132AB^. 

(2.2.13) 


The  dielectric  impermeability  matrix  now  needs  to  be  diagonalized  in 
order  to  obtain  the  major  and  the  minor  axes  of  the  cross-sectional  ellipse 
formed  by  the  plane  perpendicular  to  the  incident  light  as  it  cuts  the 
Indlcatrlx.  The  equation  of  a  quadric  surface  which  denotes  the  indicatrix 
is  given  by 


l  BkaVl  =  1  *-1.2,3) 


(2.2.14) 


where,  as  before,  Bka  =  Bkt  +  ABkt,  and  x^  x2,  and  x3  form  the  rectilinear 
coordinate  system  x,  y,  and  z  for  the  sample.  B0  is  the  isotropic  dielec- 
trie  impermeability,  given  by  BQ  =  l/n£  where  nQ  is  the  Isotropic  refractive 
Index. 

As  before,  the  problem  can  be  simplified  by  taking  the  ray  direction  to 
be  parallel  to  axis  -X3  (l.e.,  along  the  sample  -z  axis).  In  that  case, 
Intersection  of  the  plane  perpendicular  to  the  ray  direction  is  obtained  by 
setting  X3  =  0  and  Eqn.(2.2.14)  simplifies  to  give  the  cross-section  of  the 
Indlcatrlx  perpendicular  to  the  ray  direction  as 

2  2 

(Bq  +  AB11)x1  +  (Bq  +  AB22)X2  +  AB,  2xix2  +  AB^x^  =  1.  (2.2.15) 

Forming  the  matrix,  diagonalizing  It  and  solving  the  characteristic  equation, 
we  find  the  two  roots  which  are  given  by: 


AB. . +AB-o 

B1  ■  Bo +  ♦ 


i  [(4B 


11  "  AB225  +  4AB 
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12AB2lJ 


(2.2.16) 
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-  AB2Z)2  +  4AB12AB27] 


(2.2.17) 


where  Bj  and  B2  are  the  dielectric  Impermeabilities  along  the  minor  and 
the  major  axes  of  the  Index  ellipse,  respectively.  Once  they  are  known,  the 
refractive  Indices  along  those  two  directions  can  be  easily  calculated  from 
nj  =  1/VBj,  and  n£  =  I/VB2. 

The  image  Intensity  distribution  for  an  edge  dislocation  can  be  obtained 
by  following  exactly  the  same  procedure  taken  to  calculate  the  Intensity 
distribution  for  a  semiconductor  disc  under  diametrical  compression  In 
Section  2.1.  The  final  expression  for  the  Intensity  I  will  be  the  same  as 
Eqn. (2.1.12)  and  Is  rewritten  below  for  convenience. 


I  =  (A/2)  sin22(a-3)  si n2{ird(l/A1  -  1/A2) ] 


(2.2.18) 


where  A  Is  the  amplitude  of  the  Incident  wave,  a  Is  the  angle  the  polarizer 
makes  with  respect  to  the  sample  x-axIs,  p  Is  the  angle  given  by 


P  =  j  tan-1 


ABj ^  4-  AB2j 


(2.2.19) 


d  Is  the  sample  thickness,  A*  =  A/nj,  A2 
wavelength  of  the  Incident  light  In  vacuum. 


A/n2, 


It  Is  Important  to  point  out  here  that  the  approach  taken  here  to  obtain 
the  Image  Intensity  is  quite  general.  One  will  r .ed  to  obtain  the  values  for 
appropriate  coefficients  aij  in  Eqn. (2.2.6)  for  other  orientations  of  the 
edge  dislocation.  The  problem  is  considerably  more  tedious  if  the  Incident 
light  is  not  taken  to  be  along  the  sample  z  axis. 

2.3  Screw  Dislocation 

In  this  section,  an  example  of  Image  simulation  for  screw  dislocation  Is 
presented.  The  analysis  given  In  Hirth  and  Lothe  [10]  Is  followed  closely. 
Let  us  consider  a  right-handed  screw  dislocation  along  the  axis  of  a  cylinder 
of  radius  R  and  length  L,  as  shown  In  Fig. 2. 3.1.  This  dislocation  can  be 
produced  from  a  perfect  cylinder  by  shear  displacements  in  the  z  direction. 

This  dislocation  has  a  Burger's  vector  b  with  the  displacements  given  by 
uj,  U2,  and  U3  along  directions  x,  y,  and  z,  respectively.  The  latter 
directions  are  also  referred  to  as  xj,  X2,  and  X3.  The  strains  and 
stresses  are  given  by  e^j  and  o^j,  respectively,  where  both  1  and  j  take 
the  values  1,  2,  and  3.  As  Is  evident  from  this  figure,  uj  =  U2  *  0,  and 
the  displacement  U3  discontinuous  at  the  cut  surface  defined  by  y  *  0, 
x  >  0  Is  given  by 
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11m  U3(x,-Ay)  -  U3(x,Ay)  *  b  (In  2  direction);  Ay  positive 

Ay  -»  0 
x  >  0. 


Assuming  that  In  an  Isotropic  medium  the  displacement  U3  Increases 
uniformly  with  the  angle  y  to  give  the  discontinuity  represented  in  the 
cylindrical  coordinate  system  as 


Y  d  ,y  b  .  x2 

Mp.Y)  -  b  55  ■  25  tan  7  *  5  tan  ^ 


b 

2ir 


(2.3.1) 


Now  the  strains  e-j j  are  defined  by 


fc1j 

which  gives 
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(1-1. 2, 3;  J-1,2,3), 


(2.3.2) 


ell  =  e22  =  e33  =  e12  *  e21  =  °» 
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e23  5£32s  - 

^  ^  4ti  8x, 


-1  X2 
tan  1  — 


ll 


2  2  * 
4n  x.  +  x2 
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e13  = 
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-1  x2 
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2  2 
4n  xl  +  x2 


(2.3.3) 


The  stresses  and  strains  are  related  by  the  stiffness  coefficients 
c1jkH*  Conventionally  these  stiffness  coefficients  are  defined  only  along 
the  crystal  principal  axes.  However,  the  strain  elements  in 

Eqn. (2.3.2)  need  not  be  along  the  crystal  principal  axes.  The  axes  system 
chosen  to  represent  the  strains  due  to  the  screw  dislocation,  i.e. 
(x1.x2.x3)  need  not  coincide  with  any  of  the  principal  crystal  axes. 
Subsequently  In  this  section,  all  primed  notations  refer  to  the  crystal  axes 
system  and  all  unprimed  notations  refer  to  the  sample  axes  system. 

The  approach  can  be  quite  general  to  take  into  account  any  orientation 
of  the  slip  plane  with  respect  to  the  crystal  principal  axes.  As  a  specific 
example,  let  us  assume  the  sample  to  be  of  (100)  orientation,  with  the 
crystal  x‘i  axis  coming  out  from  the  sample  surface  (i.e.  along  the  sample 
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>  ffl 
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axis  X3),  and  the  crystal  x'2  axis  making  an  angle  of  iji  with  the 
dislocation  slip  plane  (l.e.  sample  X1X3  plane)  as  shown  In  Fig. 2. 3.1. 
For  this  case,  the  direction  cosines  ay  between  the  crystal  axes  and  the 
sample  axes  can  be  obtained  from  Eqns.(Z.2.6).  Once  the  direction  cosines 
are  evaluated,  the  strains  can  be  converted  from  the  sample  axes  to  the 
crystal  axes  by  the  following  transformation 


e 


kS. 


ak1alj€i j 


(1=1, 2, 3;  j=l,2,3) , 


(2.3.4) 


where  e'^j,  give  the  strain  values  In  the  crystal  axes  system.  Each  of  the 
e^  term  Is  summed  over  nine  terms.  However,  =  e22  *  €33  =  e12  ■  £2i  ■  °» 

as  well  as  e23  =  e32,  and  e31  =  e13.  Taking  these  Into  account,  Eqn. (2.3.4) 

can  be  rewritten  as 

€kH  =  (aklail3  +  ak3adl)e31  +  (ak2alL3  +  ak3aH2^e23'  (2.3.5) 

For  the  specific  case  under  consideration  for  which  the  direction  cosines  are 
given  by  Eqns. (2.2.6) ,  we  have  a, .=a19=a9,=a.„  =  0,  and  a.,=l,  and  thus 
Eqn. (2.3.5)  can  be  given  by  11  ^  1J 


ell  =  e22  =  e23  =  e32  =  e33  =  °» 

t 

e12  =  e21  =  a21e31  +  a22e23’  ! 

t 

and  • 


£13  “  e31  =  a31e31  +  a32e23*  (2.3.6) 

Now  the  stress  components  o'y  along  the  crystal  principal  axes  can  be 
evaluated  by  the  following  equation 


Oy  ■  l  c-f  j|cftekftr  (^*1.2,3;  1=1,2, 3),  (2.3.7) 

k  ,  Sr 

where  are  defined  by  Eqn. (2. 3. 6).  The  stiffness  coefficients  c^ i ^2, 

used  In  Eqn. (2.3.7)  can  be  expressed  In  reduced  notation  as  done  earlier  for 
■'Mjk!  In  Eqns. (2.1.13-2.1.15)  as 

Cy|j  =  ( 1  ♦  j t k , S=1 ,2 , 3 ;  m,n=l ,2, •  •  •  ,6) •  (2.3.8) 


As  before  for  materials  with  cubic  43m  or  m3m  symmetry,  the  following  equal¬ 
ities  are  valid: 


1 . e. ,  c 


l.e.,  c 


l.e.. 


ft 
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2  "  C43 


=  c65  "  0; 


as  there  are  only  three  non-equal  stiffness  constants  namely,  Cj.,  c12,  and 
c44<  Substituting  the  above  In  Eqn. (2.3.7)  yields 


Oii  a  *  Oi 


Oi,  *  o, 


'11  22  23  32  ■  33  "  u‘ 

’l2  *  °21  *  (c1212  +  C1221)e12  "  2c44€i2* 


°13  =  °31  =  tc1313  +  C1331)e31  *  Zc44e31* 
However,  C44  Is  the  shear  modulus  p.  Therefore, 


(2.3.9) 


"12  ■  °21  ’  2“'i2' 


and  similarly. 


°31  *  °13  *  ^31  • 

Once  the  stress  components  o'jj  are  known,  the  changes  In  the  dielectric 
Impermeability  can  be  evaluated  as  done  earlier  from 


ABA.  =  l  it 


fa  kaijvij 


(1-1. 2, 3;  j=l ,2,3) , 


(2.3.10) 


I 


where  fficllj  are  the  piezo-optic  coefficients  given  earlier  In 
Eqns. (2.1. 13-2. 1.15).  Substitution  yields. 


ABil  =  AB^2  *  “  AB32  *  AB33  “ 


AB12  ’  AB21  "  2lT1212°l2* 


and 

AB31  =  AB13  =  2lI1313°3r  (2.3.11) 

These  changes  In  the  dielectric  Impermeabll itv  are  defined  along  the 
crystal  principal  axes.  Now  they  need  to  be  transferred  to  the  sample  axes 
system  through  the  following  transformation 

AB1j  =  E  ak1aljABk!l  (k«l,2,3;  11=1,2,3).  (2.3.12) 

k ,  1 

Substitution  of  the  values  for  AB'j^  from  Eqn.(2.3.11)  and  a^  from 
Eqn. (2.2.6)  Into  Eqn.(2.3.12)  yields 

ABjj  =  AB  j  2  =  AB21  =  AB22  =  AB^^  =  0, 

AB13  =  AB31  =  a21ABl2  +  a31AB31 * 


and 

AB23  =  AB32  =  a22AB12  +  a32AB31’  (2.3.13) 

The  Indicatrix  is  given  by 

X  B..X.X.  =  1  (1-1, 2, 3;  j=l,2,3) .  (2.3.14) 

U  1  J 

Expanding  as  before  and  remembering  that  there  are  no  cross-terms  under 
unstressed  condition,  one  obtains 

(Bq  +  AB^JXj  +  (BQ  +  AB22)x2^  +  (B0  +  ab33)x32 

+  2ABj2xjX2  +  2AB23x3x3  +  2AB31x3Xj  -  1, 


or. 


2  2  2 

B„x,  +  B„x-  +  B„x,  +  2AB--X-X-  +  2AB,.x,x1  =  1. 
01  02  03  23  23  31  31 


(2.3.15) 


If  the  light  is  Incident  along  the  -X3  direction,  then  the  cross-section 
found  by  putting  X3  =  0  in  the  above  equation  gives 


.&**.**«-»• i 


»’♦  i’t  *>  h  «  *  t<|  . 


, >1  <«  t.  •  (ja  -  $««*  a*- 


and  there  Is  no  birefringence  as  the  cross-section  Is  a  circle.  Therefore, 
an  Important  observation  Is  that  screw  dislocations  cannot  be  viewed 
end-on,  l.e.,  along  the  direction  of  the  Burger's  vector.  However,  If  the 
light  traverses  along  other  directions,  one  should  be  able  to  see 
Interference  fringe  patterns  due  to  birefringence.  Assuming  the  light  to  be 
falling  along  the  -xj  direction,  the  cross-section  found  by  putting 
xj  =  0  In  Eqn. (2.3.15)  gives 


2  2 

Box2  +  Box3  +  2AB23x2x3 


(2.3.17) 


This  cross-section  Is  an  ellipse  and  one  should  see  birefringence  correspond¬ 
ing  to  the  relative  dielectric  Impermeabilities  along  the  minor  and  the  major 
axes  given  by 

B1  =  Bq  +  AB23,  and  B2  -  Bq  -  AB23 

where  AB23  Is  given  by  Eqn. (2.3.13) ,  with  corresponding  refractive  Indices 
along  Its  minor  axis  being  n1  =  l/VBj  and  along  its  major  axis  being 

n2  »  1//B2.  The  angle  p  subtended  by  the  minor  Index  axis  with  axis  x3  Is 

given  by 


P  =  -|  tan-1 


AB33  "  AB22 


(2.3.18) 


However,  AB22  a  ab33  =  0  and  Eqn. (2.3.18)  gives  the  value  of  p  as 


p  =  |  tan"1-  =  |  . 

As  before,  the  resultant  transmitted  intensity  through  the  analyzer  Is  given 
by 


I  =  (A2/2)  sin22(a-P)  s1n2{ird(l/A1  -  l/*2)} 


(2.3.19) 


where  Aj  =  A/nj,  A2  *  A/n2  and  A  Is  the  incident  light  wavelength  in 
vacuum.  Here,  a  Is  the  polarization  angle  of  the  Incident  light  measured 
with  respect  to  the  axis  X3,  and  d  Is  the  material  thickness  along  xi 
direction. 
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3.  RESULTS  AND  DISCUSSIONS 


3.1  Experimental  Set-Up 

A  dark  field  plane  polarlscope  was  constructed  to  view  the  fringe 
patterns  of  diametrically  loaded  discs  of  silicon  and  gallium  arsenide.  The 
setup  Is  shown  In  Fig. 3. 1.1.  Figure  3.1.1(a)  shows  a  general  view  of  the 
set-up  while  Fig. 3. 1.1(b)  shows  the  schematic  arrangement.  The  light  source 
used  was  a  linearly  polarized  19-mW  He-Ne  laser  (Jodon  HN-50)  tuned  to  1.15 
pi  wavelength.  Both  silicon  and  gallium  arsenide  with  absorption  cut-off 
wavelengths  of  1.11  pi  and  0.89  pi,  respectively,  are  transparent  at  1.15  pi 
wavelength.  A  Hamamatsu  infrared  camera  (model  C-1000-03)  was  used  as  a 
detector.  The  power  output  of  the  laser  was  strong  enough  for  the  detector 
to  pick  up  bright  Images.  The  laser  beam  was  spatially  filtered  and  then 
expanded  to  a  diameter  of  2.5  cm  by  a  combination  spatial  filter  -  beam 
expander  arrangement. 

A  laser  light  source  was  selected  over  the  more  conventional  filtered 
Incoherent  light  sources  because  of  Its  monochromaticity.  The  use  of  the 
laser  as  the  source  of  Illumination  proved  to  reduce  the  washing  out  effect 
observed  In  Images  formed  with  wider  bandwidth  sources.  A  comparison  of 
Figs. 3. 2. 1(a)  and  3.2.2(a)  with  Fig. 3. 2. 3(a) ,  which  Is  reproduced  here  from 
the  reported  work  of  Appel  et.al.  [11],  clearly  shows  Improved  constrast  In 
the  first  two  cases  using  a  laser  light  source.  An  unfortunate  artifact  of 
the  use  of  the  laser  Is  the  production  of  speckle  In  the  Images.  This  can  be 
seen  In  Figs. 3. 2. 1(a)  and  3.2.2(a)  as  an  overall  graininess  In  the  Images. 

In  theory  at  least  the  use  of  the  polarized  laser  beam  would  have 
precluded  the  need  for  using  a  polarizer  ahead  of  the  sample.  However,  even 
In  the  beam  expanded  state,  the  laser  source  proved  to  be  too  powerful  and 
blinded  the  Infrared  sensitive  camera.  Because  of  the  need  to  rotate  the 
plane  of  polarization  of  the  light  source,  and  also  the  need  to  reduce  the 
light  Intensity,  a  polarizer  was  Inserted  Into  the  optical  path  ahead  of  the 
sample. 

Figures  3.2.1(a)  and  3.2.2(a)  were  photographed  directly  from  the  screen 
of  the  monitor  with  a  35-mm  camera.  The  samples  were  placed  between  the 
polarizer  and  the  analyzer  and  loaded  by  a  pneumatically  operated  piston. 
The  anvils  of  the  compressor  were  made  of  flame-hardened  mild  steel  and 
pieces  of  Index  card  were  used  to  protect  the  edges  of  the  samples  during 
compression.  The  cross-section  of  the  piston  plunger  was  1  sq.  In.  ±  1%.  A 
pressure  gauge  calibrated  In  psig  then  allowed  the  compressive  force  on  the 
sample  disc  to  be  read  directly. 

3.2  Results  for  Diametrically  Compressed  Semiconductor  Discs 

The  silicon  samples  used  In  this  experiment  were  cut  along  the  (111) 
orientation  from  a  rod  of  single  crystal  silicon  of  1.746  cm  diameter.  To 
prevent  loss  of  Intensity  due  to  surface  scattering,  the  sample  surfaces  were 
polished  to  a  mirror  finish.  The  polishing  was  done  manually  on  both  sides 
of  the  samples  with  SIC  powder  on  a  glass  plate  with  light  hand  pressure, 
using  a  successively  finer  abrasive  medium  starting  with  200-grit  powder  size 
and  ending  with  1200  grit  powder  size.  After  the  mechanical  polishing  was 
completed,  the  samples  were  chemically  polished  for  2  min  In  CP-4A  solution 
containing  three  parts  HF,  five  parts  HNO3,  and  three  parts  CH3COOH  by 
volume. 


Fig. 3. 1.1(a)  The  general  view  of  the  experimental  set-up. 


I 

I 


w 

% 


SAMPLE 


■ 

I 


MONITOR 


CAMERA 


ANALYZER 


COLLIMATOR/ 

BEAM 

EXPANDER 


m 

V.' 


Ml**. 


1 


POLARIZER 


He-Ne  LASER 
(1.15  jim) 


fl 


m 


Fig. 3. 1.1(b)  The  schematic  arrangement  of  the  experimental  set-up. 
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Figure  3.2.1(a)  shows  the  observed  Image  at  1.15  pm  under  the  plane 
polarlscope  for  a  (111)  orientation  silicon  crystal  disc  cf  0.873  cm  radius 
and  0.16  cm  thickness.  The  applied  load  was  19  kg  and  the  angle  of 
polarization  a  of  the  incident  light  was  90°.  Figure  3.2.1(b)  is  the 
simulated  image  for  these  same  conditions  obtained  from  the  analysis  given  in 
Section  2.1.  The  Intensity  levels  shown  in  Fig. 3. 2. 1(b)  were  obtained  by 
normalizing  the  maximum  level  to  unity  and  then  plotting  in  sixteen  gray 
scale  levels  utilizing  a  logarithmic  scale.  Any  intensity  more  than  0.1  in 
the  normalized  scale  is  plotted  as  white  in  the  computed  figure.  Figure  3.2.2 
compares  the  observed  and  the  simulated  Images  for  the  same  conditions  as 
Fig. 3. 2.1  but  with  the  load  now  increased  to  38  kg.  A  close  observation  of 
the  simulated  image  of  Fig. 3. 2. 2  will  bring  out  the  feature  of  the  small 
lobes  around  the  point  of  loading  which  is  also  present  in  the  experimentally 
observed  image.  Thus  the  match  between  the  computer  simulations  and  the 
experimentally  observed  images  is  excellent. 

The  data  reported  by  Appel  et.al.  [11]  of  a  silicon  sample  with  radius 
R  =  0.437  cm,  thickness  d  =  0.188  cm,  and  the  polarization  angle  a  =  60°  are 
shown  in  Fig. 3. 2. 3(a) .  The  crystal  orientation  of  the  sample  was  not 
mentioned  in  their  paper.  In  this  work,  we  simulated  their  image  by  assuming 
the  sample  orientation  to  be  (100)  and  (111).  The  results  of  our  simulation 
are  shown  in  Figs. 3. 2. 3(b)  and  3.2.3(c),  respectively,  for  (100)  and  (111) 
oriented  silicon  samples.  All  other  parameters  used  in  our  simulation  were 
the  same  as  those  reported  by  Appel  et.al.  [11].  Comparing  Fig. 3. 2. 3(a)  with 
Figs. 3. 2. 3(b)  and  3.2.3(c),  an  excellent  match  is  seen  between  the  simulated 
image  and  the  experimentally  observed  image  for  a  (100)  orientation  sample, 
whereas  the  simulated  image  for  a  (111)  orientation  sample  shows  a  poor  match 
with  theexperimental  one.  The  best  match  for  the  (100)  orientation 
simulation  was  obtained  for  the  case  in  which  the  load  was  applied  along  one 
of  the  crystal  principal  axes.  Hence,  we  infer  that  the  sample  used  in  the 
work  reported  by  Appel  et.al.  [11]  was  of  (100)  orientation: 

Hence,  the  algorithm  can  not  only  be  used  to  predict  the  interference 
fringe  pattern  for  any  arbitrarily  applied  stress  direction,  but  can  also  be 
used  to  predict  the  orientation  of  an  unknown  sample  by  putting  it  under  a 
compressive  load  and  then  by  comparing  the  resulting  fringe  pattern  with  the 
simulated  ones. 

The  algorithm  was  also  used  to  determine  the  fringe  patterns  for  a 
diametrically  loaded  GaAs  disc.  The  disc  was  single  crystal  undoped 
semi-insulating  GaAs  of  (100)  orlentaton,  grown  by  the  Liquid  Encapsulated 
Czochralski  (LEC)  technique.  It  was  of  diameter  0.75  inch  and  thickness 
65  mils.  Both  its  surfaces  were  polished  to  a  mirror  finish. 
Figure  3.2.4(a)  shows  the  simulated  pattern  for  the  above  disc  under  a 
diametrical  compression  load  of  20  kg,  and  the  angle  of  polarization  a  of  the 
Incident  light  of  0#.  Figure  3.2.4(b)  shows  the  experimentally  obtained 
Image  for  the  disc  subjected  to  same  conditions  as  in  Fig. 3. 2. 4(a) .  Figures 
3.2.5(a)  and  3.2.5(b)  show  the  simulated  pattern  and  the  experimentally 
obtained  Image  of  the  same  disc  but  with  the  load  now  increased  to  40  kg  and 
the  angle  of  polarization  a  of  the  Incident  light  changed  to  45°.  The  match 
between  the  simulated  patterns  and  the  experimentally  obtained  images 
confirms  the  general  nature  of  the  algorithm  which  can  effectively  be  used 
for  a  variety  of  materials. 
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Fig.  3.2.1.  a)  Experimentally  observed  Image  on  a  (111)  silicon  disc  of 
radius  0.873  cm  and  thickness  0.16  cm.  Diametrical  compression 
load  P  »  19  kg,  angle  of  polarization  a  =  90°.  b)  Simulated 
image  for  the  same  conditions  as  a).  Normalized  intensity  Is 
plotted  on  a  logarithmic  scale  for  the  computed  image. 
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Fig.  3.2.2.  a)  Experimentally  observed  and  b)  computed  images  on  a  (111) 
silicon  disc  of  radius  0.873  cm,  thickness  0.16  cm  for  a  load  of 
38  kg  and  the  angle  of  polarization  a  of  90®.  Normalized 
Intensity  Is  plotted  on  a  logarithmic  scale  for  the  computed 
image. 
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Fig.  3.2.3.  a)  Experimentally  observed  image  reported  by  Appel  et.al.  [HJ 
for  a  silicon  sample  of  unknown  orientation  with  a  radius  of 
0.437  cm,  thickness  of  0.188  cm,  load  of  38  kg  and  the  angle  of 
polarization  a  of  the  incident  light  of  60°,  b)  simulated  image 
for  the  same  conditions  as  a)  for  (100)  orientation,  and  c)  same 
as  b)  but  for  (111)  orientation.  Normalized  intensity  is 
plotted  on  a  logarithmic  scale  for  the  computed  images. 
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Fig.  3.2.4.  a)  Simulated  Image  on  a  (100)  gallium  arsenide  disc  of  radius 
0.375  Inch  and  thickness  65  mils.  Diametrical  compression  load 
P  „  20  kg,  angle  of  polarization  a  =  0°.  b)  Experimentally 
obtained  Image  for  the  same  conditions  as  a).  Normalized 
Intensity  Is  plotted  on  a  logarithmic  scale  for  the  simulated 
Image. 
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Fig.  3.2.5.  a)  Simulated  image  on  a  (100)  gallium  arsenide  disc  of  radius 
0.375  inch  and  thickness  65  mils.  Diametrical  compression  load 
P  *  40  kg,  angle  of  polarization  a  =  45*.  b)  Experimentally 
obtained  image  for  the  same  conditions  as  a).  Normalized 
intensity  is  plotted  on  a  logarithmic  scale  for  the  simulated 
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Isotropic  materials  with  spherical  Indicatrles  become  anisotropic  when 
stressed  either  by  some  external  means  such  as  application  of  a  load  or  by 
internal  means  such  as  defects  within  the  material.  The  stress  field  causes 
optical  anisotropy  and  as  a  result  the  material  becomes  blrefrlngent  with  the 
indlcatrlx  represented  by  an  ellipsoid.  Thus,  if  the  sample  Is  Illuminated 
with  a  monochromatic  radiation  of  wavelength  longer  than  the  absorption 
cut-off  of  the  material,  upon  entering  the  stressed  sample,  the  light  will 
split  Into  two  components  along  the  major  and  the  minor  axes  of  the 

cross-sectional  ellipse  of  the  Indlcatrlx  with  different  Indices  of 
refraction.  These  two  components  will  build  up  a  phase  difference  $  between 
them  as  they  traverse  the  sample.  This  phase  difference  Is  given  by 

Cfo.-o^tid 

♦  -  - -  ,  (3.3.1) 

where  C  Is  the  stress-optic  coefficient  of  the  material,  oj,  and  02  are 

the  diagonalized  elements  of  the  stress  matrix  given  by 

„  °I!  *  °22  1  r.  „  ,2  .  ,„2  ~|l/Z  ,,  - 

®1  s  2  ‘*’2  11  ~  ®22  ^°12 J  »  (3.3.2) 

and 

qn  +  o22  1  H  ,2  .  2  “11/2  n 

°2  2  "2  L(  11  '  °22*  +  4o12  J  ’  (3.3.3) 

where  on,  022.  and  012  are  given  by  Eqns. (2. 1.1-2. 1.3)  respectively,  d 

is  the  material  thickness,  and  A  Is  the  wavelength  of  the  Incident  light. 

This  approach  assumes  a  constant  value  of  the  stress-optic  coefficient. 
However,  the  stress-optic  coefficient  can  be  taken  as  a  constant  only  for 
Isotropic  materials.  In  an  anisotropic  medium,  which  Is  the  case  for  a 
diametrically  compressed  semiconductor  sample,  the  value  of  the  stress-optic 
coefficient  C  will  depend  on  the  orientation  of  the  load  with  respect  to  the 
principal  crystal  axes.  The  values  of  C  will  also  vary,  in  general,  from  one 
point  of  the  sample  surface  to  the  next  for  a  given  load.  However,  the 

values  of  C  do  not  depend  on  the  magnitude  of  the  applied  load,  provided  that 

a  linear  relationship  between  the  pressure  and  the  phase  retardation  holds. 
Glardinl  T2]  has  reported  that  this  linearity  holds  upto  a  pressure  of 
450  kg/cm^for  silicon. 

In  this  work,  a  different  approach  is  taken  to  examine  the  behavior  of 
the  stress-optic  coefficient  C  for  the  case  of  diametrical ly  compressed  discs 
of  semiconductors  (silicon  and  gallium  arsenide).  The  Incident  light  while 
passing  through  the  stressed  sample  splits  into  two  components  with  different 
Indices  of  refraction,  and  thus  travels  with  different  velocities.  The 

velocity  of  light  within  a  material  with  Index  of  refraction  n  Is  given  by 
c/n,  where  c  Is  the  velocity  of  light  In  free  space.  Now  It  Is  assumed  that 
the  light  travelling  along  the  minor  axis  of  the  ellipse  has  a  velocity  vj, 
and  that  travelling  along  the  major  axis  has  a  velocity  V2,  with 

corresponding  wavelengths  Aj  and  A2,  and  the  Indices  of  refraction  nj 
and  n2»  The  phase  angle  of  the  minor  axis  component  Is  given  by 


♦  l  «  2nd/Ai,  and  that  of  the  major  axis  component  is  given  by 
♦2  =  2ird/A2*  The  phase  angle  difference  $  between  these  two  components 
Is  given  by 


♦  *  $2  ”  *  2ird  (^~ 


n9f  n.f 

*  =  2nd  (-§ - £-), 

*  2M  .  n  v 

*  "  A  (n2~nl* 


(3.3.4) 


where  A  Is  the  wavelength  of  the  Incident  light  In  free  space  and  f  Is  Its 
frequency. 

Comparison  of  Eqns. (3.3.1)  and  (3.3.4)  gives 


n2  "  nl 
°1  "  °2 


(3.3.5) 


Equation  (3.3.5)  was  used  to  compute  the  value  of  C  at  each  point  on  the 
wafer  for  a  semiconductor  disc  under  diametrical  compression.  For  both  SI 
and  GaAs  (100)  samples,  the  value  of  C  was  found  to  vary  along  the  sample 
surfaces.  The  value  of  C  was  also  found  to  depend  on  the  orientation  of  the 
load  with  respect  to  the  principal  crystal  axes.  However,  for  both  SI  and 
GaAs  (111)  samples,  C  was  found  to  be  Independent  of  the  loading  orientation 
or  the  position  on  the  sample  surfaces  as  expected  from  the  crystal  symmetry. 

The  values  of  C  calculated  from  the  simulation  as  a  function  of  position 
on  the  sample  surface  for  a  (100)  oriented  silicon  disc  are  shown  In 
Fig. 3. 3.1.  Figure  3.3.2  shows  the  same  for  a  (100)  oriented  gallium  arsenide 
disc.  In  both  of  these  figures,  the  orientation  of  the  load  with  respect  to 
the  principal  crystal  axes  lying  on  the  wafer  surface  is  defined  by  the  angle 
5.  Figure  3.3.1(a)  shows  the  values  of  the  stress-optic  coefficient  C  as  a 
function  of  position  on  a  (100)  SI  surface  for  S  «  45°,  while  Fig. 3. 3. 1(b) 
shows  the  same  for  S  =  75°.  Similarly,  Figs. 3. 3. 2(a)  and  3.3.2(b)  show  the 
values  of  C  for  a  (100)  GaAs  sample  for  the  angles  S  =  45°,  and  S  =  75’, 
respectively. 

From  Figs. 3. 3.1  and  3.3.2  It  can  be  seen  clearly  that,  as  expected,  C 
cannot  be  taken  as  a  constant  for  a  given  load  as  It  Is  a  function  of 
position  on  the  sample  surface.  As  the  angle  of  loading  Is  changed,  the 
pattern  of  C  on  the  wafer  changes.  Therefore,  It  Is  a  function  of  the 
orientation  of  the  loading  axes  with  respect  to  the  principal  crystal 
directions.  For  (100)  oriented  silicon,  the  values  of  C  ranged  from 
2.0  x  10“12  cm2/dyne  to  3.0  x  10“12  cm2/dyne  and  Is  plotted  In 
Fig. 3. 3.1  In  a  gray  scale  with  11  levels.  In  the  plot,  a  white  scale  was 
employed  for  C  >  2.95  x  10-12  cm2/dyne  and  a  black  scale  for 
C  <  2.05  x  10-12  cm2/dyne.  Intermediate  values  were  plotted  In  nine 
levels  varying  In  the  shade  of  gray,  with  linear  spacings.  For  (100) 
oriented  gallium  arsenide,  the  values  of  C  ranged  from  0.8  x  10-12 
cm2/dyne  to  2.6  x  10“12  cm2/dyne  and  are  plotted  In  Fig. 3. 3. 2  In  a  gray 
scale  with  10  levels.  In  the  plot,  a  white  scale  was  employed  for 
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Fig.  3.3.1  Wafer  map  of  the  computed  values  of  the  stress-optic  coefficient 
C  for  a)  6  =  45°  and  b)  5  =  75°  for  a  load  of  19  kg,  (100) 
silicon  wafer  with  radius  of  0.873  cm  and  thickness  of  0.16  cm. 
Plotted  In  11  levels  of  gray  scale  with  white  for 
C  >  2.95  x  10“12  cm2/dyne  and  black  for  C  <  2.05  x  10“12 
cm2/dyne. 
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Fig.  3.3.2  Wafer  map  of  the  computed  values  of  the  stress-optic  coefficient 
C  for  a)  6  =  45°  and  b)  S  =  75°  for  a  load  of  20  kg,  (100) 
gallium  arsenide  wafer  with  radius  of  0.375  inch  and  thickness  of 
65  mils.  Plotted  in  10  levels  of  gray  scale  with  white  for 
C  >  2.5  x  10~12  cm2/dyne  and  black  for  C  <  0.9  x  10"1Z 
cm2/dyne. 


C  >  2.5  x  10~12  cm2/dyne  and  a  black  scale  for  C  <  0.9  x  10  12  tm2/dyne.  In¬ 
termediate  values  were  plotted  in  eight  levels  varying  in  the  shade  of  gray, 
with  linear  spacings. 

One  important  point  to  note  is  that  though  C  varies  as  a  function  of 
position  for  any  given  load,  yet  for  a  given  value  of  6,  the  wafer  map  of  C 
is  independent  of  the  magnitude  of  the  applied  load.  For  instance, 
Fig. 3. 3.1  is  plotted  for  a  load  of  19  kg  and  Fig. 3. 3. 2  is  plotted  for  a  load 
of  20  kg.  However,  doubling  the  load  to  38  and  40  kg  respectively  will  not 
cause  these  two  patterns  to  change.  This  is  expected  from  the  linearity 
between  the  stress  and  the  phase  retardation. 

However,  for  (111)  oriented  Silicon  or  Gallium  Arsenide,  the  value  of  C 
is  independent  of  position  for  any  given  load,  and  is  also  independent  of  the 
orientation  of  the  loading  axes  with  respect  to  the  principal  crystal  axes. 
The  value  we  report  for  (111)  silicon  disc  under  diametrical  compression  is 
2.33  x  10-12  cm2/dyne  and  the  same  for  Gallium  Arsenide  disc  under 
diametrical  compression  is  1.94  x  10“12  cm2/dyne. 


The  value  of  C  quoted  by  Appel  et.al.  [11] 

2.1  x  10-12  cm2/dyne.  Bui  lough  r 12]  reported  it  as  variable  from 
1.0  x  10“12  cm2/dyne  to  2.3  x  10“12  cm2/dyne  based  on  his  studies  of 
edge  dislocations  in  silicon.  Prussin  and  Stevenson  [13]  deduced  a  value  of 
2.0  x  10'12  cm2/dyne  in  their  experiment  on  silicon  bars  subjected  to 
pure  bending.  Their  reported  value  compares  favorably  to  those  quoted  by 
Bui  lough  and  by  Appel  et.al.  Both  the  works  by  Appel  et.al.  and  by  Prussin 
and  Stevenson  measured  the  retardation  as  a  function  of  the  applied  stress. 
From  the  observed  linear  relation  between  them,  a  mean  value  of  the 
stress-optic  coefficient  C  was  found  to  fit  the  data,  whereas  Bui  lough 
adjusted  the  value  of  C  to  make  his  theoretical  simulation  of  edge 
dislocations  in  silicon  match  with  those  experimentally  determined  by  Bond 
and  Andrus  [14]. 

To  summarize,  we  can  say  that  for  (100)  oriented  Silicon  or  Gallium 
Arsenide  disc  under  diametrical  compression,  no  unique  value  of  C  can  be 
given  as  it  varies  from  point  to  point  for  any  given  load,  and  varies  with 
the  orientation  of  the  load  with  respect  to  the  principal  crystal  axes  as 
well.  However,  no  such  dependence  was  found  for  (111)  oriented  samples.  For 
silicon,  the  values  of  C  calculated  in  this  work  are  consistent  with  those 
reported  earlier.  However  for  gallium  arsenide  to  our  knowledge,  no  data  of 
C  is  reported  so  far.  The  value  of  C  for  a  given  position  on  the  wafer  will 
not  change  on  increasing  the  load  so  long  as  the  tensor  elements  Hijka  of 
the  piezo-optic  coefficient  matrix  remain  constant. 

3.4  Simulation  Results  for  Edqe  and  Screw  Dislocations 


for  silicon  is 
as  variable  from 
on  his  studies  of 
deduced  a  value  of 
bars  subjected  to 


The  simulation  results  for  edge  and  screw  dislocations  are  presented  in 
this  Subsection.  Figure  3.4.1  shows  the  simulation  results  for  an  edge 
dislocation  on  a  (100)  Si  surface  with  the  dislocation  slip  plane  lying  on 
one  of  the  principal  crystal  plane  and  the  Burgers  vector  along  one  of  the 
principal  crystal  direction  (i.e.,  the  angle  \|#  in  Fig. 2. 2.1  is  zero).  The 
polarization  angle  of  the  incident  light  with  respect  to  the  dislocation  slip 
plane  is  denoted  by  a.  The  images  in  Fig. 3. 4.1  are  obtained  for  the  light 
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Fig.  3.4.1  Simulated  fringe  patterns  for  an  edge  dislocation  on  a  (100)  SI 
surface  with  the  dislocation  slip  plane  along  one  of  the 

principal  crystal  plane  and  the  Burger's  vector  along  one  of  the 
principal  crystal  direction.  a  Is  the  angle  between  the 

dislocation  slip  plane  and  the  polarization  axis  of  the  incident 
light,  a)  a  =  0°,  b)  a  =  30°,  c)  a  =  45%  d)  a  =  60  . 


Incident  along  one  of  the  principal  crystal  direction.  Figures  3.4.1(a), 
(b),  (c),  and  (d)  show  the  simulation  results  for  a  =  0°,  30°,  45°,  and  60®, 
respectively.  From  the  figures,  it  can  clearly  be  seen  that  when  the 
polarizer  makes  an  angle  of  45*  with  respect  to  the  dislocation  slip  plane, 
the  image  is  brightest  with  a  four-petal  rosette  pattern  with  equal  lobes. 
This  is  of  interest  as  it  is  a  key  feature  which  distinguishes  edge 
dislocations  from  screw  dislocations  lying  along  principal  crystal  planes,  as 
will  be  evident  from  the  following  paragraphs. 

Figure  3.4.2  shows  the  simulation  results  for  a  screw  dislocation  on  a 
(100)  Si  surface  with  the  dislocation  slip  plane  lying  along  one  of  the 
principal  crystal  plane  and  the  Burgers  vector  along  one  of  the  principal 
crystal  direction.  Here  also  the  polarization  angle  of  the  incident  light 
with  respect  to  the  dislocation  slip  plane  is  denoted  by  a.  If  the  light  is 
incident  along  the  direction  of  the  Burgers  vector,  the  cross-section  is  a 
circle  and  there  is  no  birefringence  as  is  evident  from  Eqn. (2.3.16). 
Therefore,  for  the  particular  case  under  consideration,  this  dislocation 
cannot  be  viewed  end-on,  i.e.,  along  the  direction  of  the  Burgers  vector. 
Figures  3.4.2(a),  (b),  (c),  and  (d)  show  the  simulation  results  for  a  screw 
dislocation  when  the  light  is  incident  perpendicular  to  the  direction  of  the 
Burgers  vector  and  along  one  of  the  principal  crystal  direction  for  a  =  0®, 
25®,  50®,  and  75®,  respectively.  From  the  intensity  distribution  equation 
(2.3.19)  for  a  screw  dislocation,  we  can  see  that  for  a  *  p  *  45®,  there  is  a 
total  extinction  as  the  transmitted  intensity  is  zero.  However,  for  edge 
dislocations  for  a  *  45®,  the  intensity  distribution  is  brightest  with  a 
four-petal  rosette  pattern  with  equal  lobes. 

Therefore,  for  these  two  specific  cases  under  investigation,  it  is  easy 
to  find  out  the  types  of  the  dislocations,  their  locations,  and  their 
orientations  by  continuously  varying  the  polarization  angle  and  monitoring 
the  fringe  patterns  closely.  The  orientation  of  a  screw  dislocation  is  45® 
from  the  axis  of  polarization  for  total  extinction.  For  any  other  angles  of 
polarization,  the  images  of  edge  and  screw  dislocations  are  completely 
different  as  evident  upon  comparing  Figs. 3. 4.1  and  3.4.2.  However,  these 
figures  are  obtained  under  rather  restrictive  assumptions  of  the  orientations 
of  the  dislocations  and  the  incident  light  directions. 

This  simulation  approach  can  be  made  quite  general  though  it  would  be 
rather  tedious  to  work  out  for  any  arbitrary  orientation  of  the  dislocation 
slip  plane  with  respect  to  any  of  the  principal  crystal  planes,  and  any 
arbitrary  direction  of  the  Burgers  vector  with  respect  to  the  dislocation 
slip  plane.  The  approach  would  include  determining  the  direction  cosines 
between  the  dislocation  slip  plane  and  the  principal  crystal  directions. 
Then  any  arbitrary  orientation  of  the  Burgers  vector  with  respect  to  the 
dislocation  slip  plane  can  be  taken  care  of  by  resolving  it  along  the 
dislocation  slip  plane  and  perpendicular  to  It  which  would  then  give  the  edge 
and  the  screw  components.  Then  the  stress  components  form  each  of  these  two 
types  of  defects  need  to  be  added  keeping  in  perspective  the  directions  of 
the  principal  stress  components  and  the  shear  stress  components.  Then  these 
added  stress  components  need  to  be  transformed  to  the  principal  crystal 
directions  from  the  known  values  of  the  direction  cosines  between  the 
principal  crystal  directions  and  the  dislocation  slip  plane.  After  that,  the 
changes  in  the  dielectric  impermeability  can  be  obtained  from  the  known 
values  of  the  piezo-optic  coefficients. 


Fig. 3. 4. 2  Simulated  fringe  patterns  for  a  screw  dislocation  on  a  (100)  SI 
surface  with  the  dislocation  slip  plane  along  one  of  the 

principal  crystal  plane  and  the  Burger's  vector  along  one  of  the 
principal  crystal  direction.  a  Is  the  angle  between  the 

dislocation  slip  plane  and  the  polarization  axis  of  the  Incident 
light  a)  a  =  0‘,  b)  a  =  25%  c)  a  =  50%  d)  a  =  75°. 


Then  we  need  to  fix  the  direction  of  the  Incident  light  with  respect  to 
these  principal  crystal  directions  and  obtain  the  direction  cosines  between 
them.  After  that  the  changes  In  the  dielectric  Impermeability  obtained  In 
the  principal  crystal  axes  need  to  be  transformed  to  the  coordinate  system 
determined  by  the  light  direction  from  the  known  values  of  the  direction  cos¬ 
ines.  The  diagonal izat Ion  of  the  dielectric  Impermeability  matrix  then 
follows  to  obtain  the  refractive  Indices  along  the  major  and  the  minor  axes 
directions  and  then  the  procedure  of  obtaining  the  Intensity  of  the  light 
transmitted  through  the  analyzer  Is  the  same  as  outlined  In  either  of  the 
Sections  2.2,  or  2.3. 


Thus,  the  generalization  of  the  problem  is  straightforward  In  concept 
but  the  algebra  Involved  Is  sufficiently  complicated.  Currently  work  Is 
being  done  to  make  this  algorithm  more  general  In  approach.  A  variety  of 
these  Images  for  varying  directions  of  the  Burgers  vector  with  respect  to  the 
dislocation  slip  plane  and  varying  orientations  of  the  dislocation  slip  plane 
with  respect  to  the  principal  crystal  planes  for  varying  directions  of  the 
Incident  light  and  varying  angles  of  polarization  can  be  generated  and  stored 
on  the  computer.  Then  a  sample  with  any  arbitrary  direction  of  the  Burgers 
vector  and  any  arbitrary  orientation  of  the  dislocation  slip  plane  can  be 
placed  under  the  polarlscope  and  the  experimental  image  obtained  for 
different  directions  of  the  Incident  light  and  different  angles  of 
polarization.  Then  an  Image  matching  procedure  Is  needed  between  the 
simulated  Images  and  the  experimentally  observed  ones  to  determine  the  exact 
direction  of  the  Burgers  vector  and  the  exact  orientation  of  the  dislocation 
slip  plane. 
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4.  SUMMARY  AND  CONCLUSIONS 
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A  fast,  accurate,  and  efficient  technique  which  can  characterize  large 
defects  In  each  semiconductor  substrate  prior  to  device  fabrication  In  a 
reasonable  time  will  Improve  the  yield  of  the  electronic  devices  fabricated 
on  It.  In  this  work,  the  technique  of  Infrared  piezo-birefringence  is 
Investigated  to  characterize  large  defects  in  semiconductor  materials.  A 
dark  field  plane  polarlscope  using  1.15  pm  wavelength  He-Ne  Infrared  laser  as 
the  light  source  was  constructed  for  this  purpose. 

As  a  first  step  towards  defect  characterization,  the  cases  of  diametri¬ 
cally  loaded  discs  of  semiconductor  samples  were  taken.  An  algorithm  was 
developed  to  be  used  In  the  computer  for  the  simulation  of  the  stress 
patterns  and  the  resulting  fringe  patterns  for  Infrared  transparent  semicon¬ 
ductor  discs  under  diametrical  compression  placed  In  a  dark  field  plane 
polarlscope  between  crossed  polarizer  and  analyzer.  This  was  done  to  obtain 
a  better  understanding  of  the  simulation  algorithm  and  subsequently  for  Its 
use  for  defect  characterization.  The  simulation  results  matched  remarkably 
well  with  the  experimentally  obtained  Images  for  diametrically  compressed 
discs  of  SI  and  GaAs. 

The  behavior  of  the  stress-optic  coefficient  C  was  also  Investigated  In 
this  work,  which  had  been  taken  as  a  constant  In  earlier  works  by  other 
Investigators.  Our  Investigation  showed  that  for  (100)  oriented  SI  and  GaAs 
discs  under  diametrical  compression,  the  value  of  C  changed  from  one  point  to 
the  next  on  the  sample  surface  for  a  given  load.  The  value  of  C  was  also 
found  to  depend  on  the  orientation  of  the  applied  load  with  respect  to  the 
principal  crystal  axis.  However,  the  value  of  C  on  any  posltl  •’  on  the  wafer 
did  not  change  upon  Increasing  the  magnitude  of  the  load.  This  Is  expected 
from  the  linearity  between  the  stress  and  the  phase  retardation.  The  values 
of  C  for  both  SI  and  GaAs  discs  of  (111)  orientation  were  found  to  be  inde¬ 
pendent  of  position  for  any  given  load,  as  well  as  were  Independent  of  the 
orientation  of  the  load  with  respect  to  the  principal  crystal  axis,  as 
expected  from  the  crystal  symmetry.  The  values  we  report  for  (100)  SI  disc 

-12  2 

subjected  to  diametrical  compression  ranged  from  2.0  x  10  cm  /dyne  to 
-12  2 

3.0  x  10  cm  /dyne,  and  for  (100)  GaAs  disc  subjected  to  diametrical 
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compression  ranged  from  0.8  x  10  cm  /dyne  to  2.6  x  10  cm  /dyne.  For 
(111)  oriented  SI  and  GaAs  discs  under  diametrical  compression,  the  values 

are  2.33  x  10-12  cm2/dyne  and  1.94  x  10~12  cm2/dyne,  respectively. 

Next  the  algorithm  was  used  to  determine  the  fringe  patterns  for  edge 
and  screw  dislocations.  For  simplification  purposes.  Image  simulations  were 
carried  out  for  edge  and  screw  dislocations  when  the  dislocation  line,  the 
Burgers  vector  and  the  direction  of  the  Incident  light  were  along  the  crystal 
principal  axes.  For  this  special  case  when  the  Incident  light  was 
perpendicular  to  the  Burgers  vector,  the  image  patterns  for  an  edge  and  a 
screw  dislocation  were  obtained  and  they  showed  a  distinguishing  difference 
between  them.  For  the  polarization  angle  of  the  light  of  45*  with  respect  to 
the  dislocation  slip  plane,  the  Image  of  a  screw  dislocation  was  found  to 
disappear,  whereas  for  an  edge  dislocation  for  the  same  polarization  angle, 
the  Image  pattern  was  the  brightest  with  a  four-petal  rosette  pattern  with 
equal  lobes.  This  feature  was  striking  as  It  gave  a  distinguishing  factor 
for  detecting  a  screw  dislocation  from  an  edge  dislocation  by  continuously 
varying  the  polarization  angle  and  monitoring  the  fringe  patterns. 
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Though  specific  cases  were  Investigated  here  for  the  edge  and  the  screw 
dislocations,  the  approach  can  be  made  quite  general  to  take  Into  account  any 
arbitrary  direction  of  the  Burgers  vector,  the  dislocation  slip  plane,  the 
direction  of  the  Incident  light,  and  the  polarization  angle  of  the  Incident 
light.  The  concept  Is  rather  straightforward  but  the  algebra  involved  is 
sufficiently  complicated  and  tedious.  The  general  appraoch  Is  given  below. 

First,  the  direction  cosines  between  the  dislocation  slip  plane  and  the 
principal  crystal  directions  are  obtained.  Then  any  arbitrary  direction  of 
the  Burgers  vector  with  respect  to  the  dislocation  slip  plane  Is  taken  Into 
account  by  resolving  it  along  the  dislocation  line  and  perpendicular  to  It, 
which  are  the  screw  and  the  edge  components.  Then  the  stress  components  for 
these  two  cases  are  added  keeping  In  mind  the  actual  directions  of  the 
principal  stress  components  and  the  shear  stress  components.  After  that  the 
added  stress  components  are  transformed  to  the  crystal  principal  directions 
from  the  known  values  of  the  direction  cosines  between  the  principal  crystal 
directions  and  the  dislocation  slip  plane. 

Then  the  changes  in  the  dielectric  impermeability  along  the  principal 
crystal  directions  are  obtained  from  the  known  values  of  the  piezo-optic 
coefficients  and  the  stress  components.  Then  the  direction  of  the  Incident 
light  needs  to  be  fixed  with  respect  to  the  principal  crystal  directions  and 
obtain  the  direction  cosines  between  them.  After  that,  the  changes  In  the 
dielectric  Impermeability  is  transformed  to  the  coordinate  system  defined  by 
the  light  direction  from  the  known  values  of  the  direction  cosines.  Then  the 
diagonal Izatlon  of  the  dielectric  Impermeability  matrix  follows  to  obtain  the 
refractive  Indices  along  the  major  and  the  minor  axes  directions.  After 
that,  the  usual  procedure  of  obtaining  the  Intensity  of  the  transmitted  light 
through  the  analyzer  as  outlined  In  Sections  2.2,  or  2.3,  can  be  utilized. 

Currently  work  Is  being  done  to  generate  the  Images  of  edge  and  screw 
dislocations  for  varying  directions  of  the  Burgers  vector,  the  dislocation 
slip  plane,  the  Incident  light  direction,  and  the  polarization  angle  of  the 
Incident  light.  The  simulation  results  can  then  be  stored  in  a  computer. 
Then  a  sample  with  any  arbitrary  direction  of  the  Burgers  vector  and  any 
arbitrary  orientation  of  the  dislocation  slip  plane  can  be  placed  under  the 
polarlscope  and  the  experimental  Images  obtained  for  varying  directions  and 
angles  of  polarization  of  the  Incident  light.  Then  an  image  matching 
procedure  can  be  utilized  between  the  simulated  images  and  the  experimentally 
observed  ones  to  determine  the  exact  direction  of  the  Burgers  vector  and  the 
exact  orientation  of  the  dislocation  slip  plane. 
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6.  APPENDICES 


6.1  Direction  Cosines  for  a  (100)  Plane  under  Arbitrary  Loading 


Figure  6.1.1  shows  the  axes  system  chosen.  The  loading  axes  system  is 
(Xj.X2.X2)  and  the  crystal  principal  axes  system  Is  (x^.x^.xp.  Axis  X3  In 

the  loading  axes  system  and  axis  x'j  In  the  crystal  principal  axes  system 
come  out  of  the  plane  of  the  paper.  The  crystal  principal  axis  x'2  makes 
an  angle  a  with  axis  of  the  loading  axes  system.  The  load  P  Is  applied 
along  axis  X£.  The  coordinates  being  so  defined,  the  direction  cosines  can 
now  be  obtained.  These  are  given  by 


a12 

a13 

a21 

a22 

a23 

a31 

a32 

a33 


COS  Xj^Xj  = 
COS  Xj^  = 
COS  Xj^X-j  = 
COS  X^Xj  = 
COS  x£~X2  = 
COS  X^X^  = 
COS  X^Xj  = 
COS  X^X2  = 
COS  X^X3  = 


cos  90  =  0, 
cos  90  =  0, 
cos  0=1, 
cos  a, 
cos  (90-a), 
cos  90  =  0, 
cos  (90+a), 
cos  a, 
cos  90  =  0. 


(6.1.1) 


6.2  Direction  Cosines  for  a  (111)  Plane  under  Arbitrary  Loading 


The  case  for  (111)  plane  Is  more  involved.  Here,  (x1.x2.x3) 
coordinate  system  refers  to  the  loading  axes  system  with  X2  being  the  axis 
of  loading.  Axis  xi  lies  on  the  sample  surface  while  axis  X3  Is  normal 
to  the  plane  of  the  paper  as  shown  In  Fig. 6. 2.1.  Also,  (xj.x^.x^)  refers  to 

the  crystal  principal  axes  coordinate  system  as  shown  In  Fig. 6. 2. 2.  The 
origin  of  both  these  coordinate  systems  Is  taken  to  be  the  substrate  center 
0. 


Let  (x1.x2.x3)  be  the  point  on  the  wafer  periphery  where  the  load 
Is  applied.  This  arbitrary  loading  direction  can  be  specified  by  an  angle  y 
measured  in  the  clockwise  direction  from  axis  x"2  as  shown  In 
Fig. 6. 2. 3(a).  The  latter  Is  the  projection  of  crystal  principal  axis  x'i 
on  the  wafer  plane  as  Indicated  In  Fig. 6. 2. 3(b) .  Axes  X3  and  x"3  are  the 
same  axis  normal  to  the  plane  of  the  paper.  Axis  x"i  is  In  the  plane  of 
the  wafer  and  Is  perpendicular  to  both  x"2  and  x"3.  We  need  to  determine 
the  direction  cosines  a-j  j  that  relate  a  point  (x1.x2.x3)  In  the 
loading 
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axes  system  to  a  point  (x|,x^,xp  In  the  crystal  axes  system.  The  transfor¬ 
mation  equation  Is  given  as 


=  Z  ^1j*j  (J=l»2»3). 

This  linear  transformation  Is  carried  out  In  two  steps  as  given  below 


(6.2.1) 


XC  =  £  skjxj 


xi  =  £  ai’kxk 
xi  *  l  afk<I  akjxJ>' 


Comparing  Eqns. (6.2.1)  and  (6.2.2)  gives 


(J-l.2,3), 


(k=l,2,3) , 


au  ■  l  ai'kakj 


(k-1,2,3). 


(6.2.2) 


(6.2.3) 


From  Fig. 6. 2. 3,  one  obtains  by  direct  observation  that 
a^  *  cos(xpXj)  =  cos  y. 

aj2  *  cos(xpx2)  =  cos(90°-y), 
a'u  *  cos(x^x3)  *  0, 

&2i  *  costx^Xj)  =  cos(90#+y), 
a^2  =  cos(x^x2)  =  cos  Yt 
a23  =  cos(xJpx3)  =  0, 

*3!  =  cos(x^x1)  =  0, 
a^2  =  cos(x^x2)  =  0, 
a^3  =  cos(x^x3)  =  1. 


(6.2.4) 


In  the  crystal  axes  system,  from  Fig. 6. 2. 2,  we  can  write  the  unit  vector 
x2"  along  the  axis  xjj 


A  O  A  1  >N  1  ^ 

y  II  .  -  Y  I  ±  v  I  __  X-  y  I 

2  "  J~e  1  Vs  2  ^6  3 


(6.2.5) 
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where  Xj',  x2‘,  x3'  are  unit  vectors  along  the  three  crystal  principal  axes. 
This  Is  obtained  by  considering  the  Intercept  of  the  (111)  crystal  plane 
on  x^  axis  as  (a, 0,0)  and  then  by  translating  the  origin  by  (|,  |,  |)  to 

the  substrate  center  to  give  the  point  coordinates  as  (^|,  -  -  |).  Here, 
a  Is  the  crystal  lattice  constant.  The  unit  vector  pointing  from  the 
origin  (0,0,0)  to  this  point  Is  given  by  Eqn.(6.2.5).  The  unit  vector  x3" 

along  the  direction  perpendicular  to  the  (111)  plane  Is  given  from  analytical 
geometry  as 


A 


(6.2.6) 


The  unit  vector  Xj"  =  x2"  *  x3"  and  Is  given  by 


A 


(6.2.7) 


Using  Eqns. (6. 2. 5-6. 2. 7)  gives 
ajj  =  cos(x^xp  *  0, 


aj2  -  cos(x’"xp  -  cos(35.264°)  *  2/76, 
aj3  =  cos(xpc§)  =  cos(54.736°)  =  1/73, 
a^  =  cos(x^xj)  *  cos(135°)  =  -  1/72, 
a^2  =  cos(x^x")  =  cos(114.095e)  =  -  1/76, 
a23  *  cos(x£~xp  =  cos(54.736#)  =  1/73, 
a31  =  cos(x^xj)  =  cos(45°)  =  1/72, 
a32  -  cos(x£x£)  =  cos(114.095#)  =  -  1/76, 
a^3  =  cos(x^xp  *  cos(54.736*)  =  1/73. 


(6.2.8) 


Using  Eqns. (6.2.4)  and  (6.2.8),  we  can  get  the  direction  cosines  ai j 
from  Eqn .(6.2.3).  These  are  given  by 
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(6.2.9) 


For  the  special  case  when  the  loading  axis  x2  coincides  with  the  axis 
x2,  y  =  0  and 

all  =  °*  a12  =  j=^  ’  a13  =  ’  a21  =  "  ’  a22  =  "  * 

a23  '  ^  •  a31  =  /=  •  a32  =  -  ^  •  amJ  *33  ‘  /j  • 

For  any  other  value  of  y.  the  direction  cosines  are  given  by  Eqn. (6.2.9). 

6.3  Plezo-ODtlc  Coefficients  for  Gallium  Arsenide  from  the  Elasto-Optlc  and 


the  Stiffness  Coefficients 


In  this  section,  the  piezo-optic  coefficients  for  GaAs  are  derived  from 
the  elasto-optic  and  the  stiffness  coefficients.  The  notations  used  by 
Nye  [8]  are  followed.  The  compliance  Is  denoted  by  s,  the  stiffness  by  c, 
the  elasto-optic  coefficients  by  p,  and  the  piezo-optic  coefficients  by  it. 
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All  of  these  parameters  are  fourth-rank  tensors  with  each  of  them  having  a 
total  of  possible  81  terms.  The  generalized  form  of  Hooke's  law  can  be 
written  as: 

€1j  =  E  S1  jki0k2,  (k=l,2,3;  #.=1,2,3) ,  (6.3.1) 

k  ,  It 

where  k  and  8,  are  summation  Indices,  the  eij  denotes  strain,  the 
denotes  stress,  and  the  denotes  crystal  compliance.  Equation  (6.3.1) 

can  also  be  written  In  an  alternate  form  as 

°1j  =  E  C1  jk#.ek!t  (k=1*2*3;  8-al  ,2,3) ,  (6.3.2) 

k ,  & 

where  c^j^  are  the  stiffness  constants  of  the  crystal.  For  both 
and  Cijicj,,  the  following  relations  are  true 

S1  jkl  =  si j#.k  “  sjikJL  =  sj19.k* 
and 


ci jkt  =  ci jtk  a  cj1kSL  *  cji#.k* 


(6.3.3) 


Using  the  reduced  notation  of  Nye  as  given  In  Table  2.1.1,  Eqns. (6.3.1)  and 
(6.3.2)  can  be  written  as 


et  *  z  sij°j 


(1 , j=l,2,. . . ,6) , 


and 


°i  =  X 


(i , j=l ,2, ... ,6) . 


(6.3.4) 


However,  changing  the  notation  from  tensor  to  matrix  introduces  factors  of  2 
and  4  as  follows.  For  the  compliances: 


s1jkl  =  smn  when  m  and  n  are  2  or  3* 

2si jk#,  *  smn  wf,en  e1ther  m  or  n  are  5  or 

4si jkB.  “  smn  when  m  and  n  are  ®  or  (6.3.5) 

However,  for  the  stiffness  no  factors  of  2  or  4  are  necessary.  It  can  simply 
be  written  as: 


c1jkH  *  cmn  (1 ,j,k,i=l,2,3;  m,n=l,2,... ,6). 


(6.3.6) 


The  total  maximum  possible  number  of  non-equal  elements  after  taking  Into 
account  the  Identities  (6.3.3)  Is  36.  However,  as  the  crystal  possesses  more 
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and  more  symmetry  elements,  the  number  of  non-equal  parameters  reduces.  For 
the  cubic  crystals,  this  number  becomes  only  three  with  the  matrix  elements 
given  by 
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(6.3.7) 


where  x  can  stand  for  c,  s,  p,  or  m. 


For  cubic  systems,  s^j,  S12,  S44  and  cn,  cj2»  C44,  are 

related  by  [8] 


S11  +  s12 


(S1  (S1  1+2S1?)  H* 


’11  *i2M*irMi2' 

_ 1 !i2 _ c 

(Sn  -  s12)(s11+2s12)  -  C12* 

1 


44  "  c44  ‘ 

Solving  the  above  equations  for  sn  and  S12,  we  get: 

«  C11  I  C12 

n  =  (CU"C12)(C11+2C12)  ’ 

s  =  _ lll2 _  # 

12  (CU"C12)(C11+2C12) 


(6.3.8) 


(6.3.9) 


In  the  absence  of  any  applied  electric  field,  the  change  In  the 
dielectric  Impermeability  due  to  the  strain  can  be  given  as  [8] 


B  Z  pi jrs€rs  (r=l ,2,3;  s=l,2,3), 
r ,  s 


where 


PlJrs  =k^  ^ijk^kHrs  (I(b1,2,3;  *-=1.2,3), 


W1  jkH  a  E  P1  j rssrskJL  (r=1,2,3;  s=l,2,3). 

1*5 


(6.3.10) 


(6.3.11) 


(6.3.12) 
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The  elasto-optlc  coefficients  pyrs  are  dimensionless,  whereas  the  piezo¬ 
optic  coefficients  have  the  dimension  of  stress.  Both  it  and  p  follow 

the  same  matrix  form  in  reduced  notation  as  given  In  Eqn. (6.3.7). 

Equation  (6.3.12)  can  be  written  In  the  reduced  form  as: 

it  ■  T  p  s  (r*l,2,...,6),  (6.3.13) 

mn  £  Hmr  rn  1  ’  .  '*  v  ' 

where  m  and  n  can  take  on  values  between  1  and  6.  However,  for  cubic 
systems,  as  Is  the  case  here,  the  only  possible  values  of  non-zero  and 
non-equal  components  are  iru,  ni2»  and  *44*  The  non-zero  components 
are  defined  by  the  Identities 

1)  iTjj  =  ir22  =  tt33, 

2)  tt12  *  tt23  =  tt31  =  tt13  =  tt32  «  it21, 

3)  *44  =  *55  =  *66* 

The  rest  of  the  components  are  zero. 

From  Eqn. (6.3.13) 

X  Mlrsrl  (r=l ,2, . . . ,6) 


which  reduces  to 


*44  =  p44s44* 


(6.3.16) 


The  values  of  p^,  p^.  and  P44,  as  well  as  of  c^,  C12.  and 
C44  are  given  In  the  paper  by  Booyens  and  Basson  [9]  and  are  shown  In 
Table  6.3.1. 

Table  6.3.1 

The  values  of  the  elasto-optlc  coefficients  (p^j)  and  the 
stiffness  (c-ji)  for  GaAs  after  Booyens  and  Basson  [9] 

(Pi j  are  dimensionless,  c^j  are  In  dyne/cm2) 


cn 

c12 

C44 

118.8  x  1010 

53.8  x  1010 

59.4  x  1010 

The  coefficients  c^j  are  easily  converted  to  the  sjj  by  Eqns. (6.3.8) 
and  (6.3.9).  This  yields  sn  =  1.1728  x  10-12  cm2/dyne,  S12  = 
-3.6558  x  10"12  cm2/dyne,  and  S44  *  1.6835  x  10~12  cm2/dyne. 

Substituting  for  the  values  of  sy  obtained  above  and  the  values  of  p^j 
from  Table  6.3.1  in  Eqns. (6.3. 14-6.3. 16) ,  one  obtains  the  piezo-optic 
coefficients  as  follows: 

*  -  9.115  x  10“14  cm2/ dyne, 

1 A  O 

tt19  =  -  5.269  x  10”1*  cnr/dyne. 


(6.3.17) 


ir44  =  -  1.212  x  10“ 13  cm2/dyne. 

However,  for  the  TTij,  the  following  Identities  hold 
*mn  ’  ”ljkl  k,hen  n-1-2-3- 


’mn  ’  2"ljkl  whe"  "*4-5-6- 

Taking  Into  account  the  above  Identities,  the  piezo-optic  coefficients  for 
GaAs  are  given  as 

-14  2 

*1111  =  “  9. 115  x  10  cm  /dyne, 

-14  2 

1I1122  =  “  5,269  x  10  cm  /dyne, 


-14  2 

^1212  =  '  6,061  x  10  cm  /dyne. 


(6.3.18) 
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